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ABSTRACT 

Goddard High Resolution Spectrograph (GHRS) observations have revealed the presence of broad wings in 
the transition-region lines of AU Mic and Capella. It has been proposed that these wings are signatures of 
microflares in the transition regions of these stars and that the solar analog for this phenomenon might be the 
“transition region explosive events” discussed by Dere, Bartoe, & Brueckner. We have analyzed GHRS 
observations of Procyon (F5 IV-V) and HR 1099 (K1 IV + G5 IV) to search for broad wings in the UV 
emission lines of these stars. We find that the transition-region lines of HR 1099, which are emitted almost 
entirely by the K1 star, do indeed have broad wings that are even more prominent than those of AU Mic and 
Capella. This is consistent with the association of the broad wings with microflaring since HR 1099 is a very 
active binary system. In contrast, the transition-region lines of Procyon, a relatively inactive star, do not show 
evidence for broad wings, with the possible exception of N v a 1239. However, Procyon’s lines do appear to 
have excess emission in their blue wings. 

Linsky et al. found no evidence for broad wings in Capella’s chromospheric lines, but we find that the Mg n 
resonance lines of HR 1099 do have broad wings. The striking resemblance between HR I099’s Mg n and 
C iv lines suggests that the Mg n line profiles may be regulated by turbulent processes similar to those that 
control the transition-region line profiles. If this is the case, microflaring may be occurring in the Kl star’s 
chromosphere as well as in its transition region. However, radiative transfer calculations suggest that the 
broad wings of the Mg n lines can also result from normal chromospheric opacity effects rather than pure 
turbulence. The prominence of broad wings in the transition region and perhaps even chromospheric lines of 
active stars suggests that microflaring is very prevalent in the outer atmospheres of active stars. 

Subject headings: line: formation — stars: chromospheres — stars: individual (a Canis Minoris, HR 1099) 
ultraviolet: stars 


1. INTRODUCTION 

The Goddard High Resolution Spectrograph (GHRS) 
aboard the Hubble Space Telescope ( HST ) provides the sensi- 
tivity and spectral resolution necessary to study UV emission- 
line profiles that originate from stellar chromospheres and 
transition regions in much greater detail than heretofore pos- 
sible. For example, using the GHRS, Linsky & Wood (1994) 
discovered that the C iv and Si iv lines of the very active dMOe 
star AU Mic are not Gaussian in shape but possess very broad 
wings. They hypothesized that the broad wings are signatures 
of microflares in the transition region of AU Mic and proposed 
that the solar analog for this phenomenon might be the 
“transition region explosive events” discussed by Dere, 
Bartoe, & Brueckner (1989). Linsky & Wood estimated the line 
flux originating from the microflares by fitting two Gaussians 
to the emission lines— a broad component representing the 

1 Based on observations with the NASA/ESA Hubble Space Telescope , 
obtained at the Space Telescope Science Institute, which is operated by the 
Association of Universities for Research in Astronomy, Inc., under NASA 
contract NAS5-26555. 

2 Staff member. Quantum Physics Division, National Institute of Standards 
and Technology, Boulder, CO 80309-0440. 


microflaring contribution and a narrow component represent- 
ing the quiescent plasma’s emission. The broad component 
accounted for ^40% of the C iv and Si iv line flux for AU Mic, 
whereas transition-region explosive events produce no more 
than 5% of the Sun’s C iv and Si iv line flux (Dere et al. 1989). 
Compared to the Sun, the broad-component emission on AU 
Mic accounts for a far larger percentage of the transition- 
region line flux, which is reasonable because AU Mic is much 
more active than the Sun. 

Broad wings have also been found in the transition-region 
lines of Capella (Linsky et al. 1995b). Because the G1 III star 
that dominates Capella’s UV line emission is a very active star, 
Capella’s transition-region line profiles were interpreted and 
modeled in the same manner as AU Mic’s lines. The broad 
component’s contribution to the line flux is correlated with 
line formation temperature, with this contribution ranging 
from 52% for N v 21239 (log T = 5.25) to 24% for Si hi 21207 
(log T = 4.7). No broad wings were detected in chromospheric 
lines such as the O i 21300 triplet (log T ~ 4.0). Linsky et al. 
also noted that separating the Gl star’s emission into broad 
and narrow components (and subtracting the emission from 
Capella’s G8 III star) results in a smoother emission-measure 
distribution for the narrow component of the Gl star. The 
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emission-measure distribution of the broad component was 
found to resemble that of a sunspot. 

In this paper we present and analyze GHRS observations of 
Procyon and HR 1099. One of our goals is to search for broad 
components in the UV emission lines of these stars. The 
hypothesis that the broad components are associated with 
microflaring suggests that the broad line wings are a manifesta- 
tion of stellar activity and are prominent in the spectra of AU 
Mic and Capella because these are active stars. We would not 
expect to see broad wings in Procyon’s transition-region lines 
that are as prominent as those seen for AU Mic and Capella 
because Procyon is significantly less active than these stars and 
is only a little more active than the Sun (Ayres, Marstad, & 
Linsky 1981). Therefore, our Procyon data will test the hypoth- 
esis that broad transition-region line wings are caused by 
microflaring. 

Procyon ( = a CMi = HD 61421) is a nearby (d — 3.5 pc) F5 
IV-V star with a mass of 1.75 M© and a radius of 2.1 R Q 
(Strand 1951; Steffen 1985; Irwin et al. 1992). The photospheric 
properties of Procyon have been exhaustively studied (e.g., 
Griffin 1971 ; Steffen 1985; Nordlund & Dravins 1990; Dravins 
& Nordlund 1990a, b), and Brown & Jordan (1981) have used 
IUE data to construct models of Procyon’s chromosphere and 
transition region. An extensive, 24 day IUE observing cam- 
paign failed to detect any significant UV line flux variations, 
which suggests that the regions responsible for Procyon’s UV 
emission lines uniformly cover the stellar surface (Ayres 1991). 
The coronal emission from Procyon has been measured by 
many X-ray and EUV instruments, including those aboard 
Einstein , EX OS AT, ROSAT , and the Extreme-Ultraviolet 
Explorer (EUV E )( Lemen et al. 1989; Wood et al. 1994; Drake, 
Laming, & Widing 1995). 

HR 1099 ( = V71 1 Tau = HD 22468) is an RS CVn binary 
system whose components have spectral types K1 IV and G5 
IV, masses of 1.4 and 1.1 M 0 , and radii of 3.9 and 1.3 R Q , 
respectively. This star system, which is located 36 pc away, has 
an orbital period of 2.84 days (Fekel 1983; Strassmeier et al. 
1988). Optical spectroscopy and Doppler imaging techniques 
have revealed the presence of dark spots on the surface of the 
K1 star, which have been found to evolve with time (Bartolini 
et al. 1983). Ayres & Linsky (1982) found that the K.1 star is 
responsible for most of HR 1099’s UV emission-line flux. 
Further IUE observations have revealed variability in those 
lines on timescales ranging from days to years (Dorren & 
Guinan 1990). Analysis of HR 1099’s extremely strong X-ray 
emission has revealed the presence of coronal temperatures at 
least as high as log T — 7.3 (Schmitt et al. 1990; Dempsey et al. 
1993; Brown 1994). 

Because of their rapid rotation, induced by tidal locking, RS 
CVn systems as a group are extremely active in terms of bright 
UV emission lines and variability, including flaring. HR 1099 is 
one of the most active of the RS CVn binaries and therefore 
represents an extreme of stellar activity. If stellar activity does 
indeed produce “broad components’ 1 in transition-region 
lines, one would anticipate that the broad components of HR 
1099’s lines should be at least as prominent as those found for 
AU Mic and Capella. Like Procyon, HR 1099 will therefore 
provide a test of the activity-broad component correlation. 
The desire to contrast the line profiles of a relatively inactive 
star with those of a very active star is the primary motivation 
for discussing the Procyon and HR 1099 observations in the 
same paper. In addition to the search for broad line wings, we 
will also discuss other interesting aspects of these data, includ- 
ing line redshifts and density diagnostics. 


2. GHRS OBSERVATIONS OF PROCYON AND HR 1099 

The GHRS observations of Procyon and HR 1099 are sum- 
marized in Table 1. (See Brandt et al. 1994 and Heapet al. 1995 
for a description of the GHRS instrument and its capabilities 
prior to the installation of COSTAR.) In addition to the data 
listed in Table 1, several other GHRS observations of HR 1099 
were made during the week of 1993 September 14-21, which 
provide more orbital phase coverage and provide the data 
necessary to investigate line variability. These additional data 
are discussed in detail by Dempsey, Neff, & Linsky (1996). All 
the observations of Procyon and HR 1099 were made prior to 
the 1993 December installation of COSTAR into HST , which 
significantly improved the spectral resolution of GHRS, espe- 
cially for observations made through the large science aperture 
(LSA). The approximate spectral resolution of each of our pre- 
COSTAR observations is listed in Table 1. 

Images of the platinum spectral-calibration lamp were taken 
in the wavelength regions of interest prior to most of our 
spectra to maximize the accuracy of our wavelength cali- 
brations. The exceptions are the HR 1099 observations of the 
1637-1672 and 1532-1568 A spectral regions. For these 
spectra, the standard wavelength calibrations were used, but 
with offset corrections based on a wavelength-calibration 
image made before the source spectra with the same grating 
but of a different spectral region — a “ spectrum y-balance,” or 
SPYBAL, calibration image. (SPYBAL calibration spectra are 
used to properly center spectra on the diode array and are 
usually made whenever the observed wavelength region is 
changed.) The dispersion relations of the SPYBAL-calibrated 
spectra are roughly 3-10 times less accurate than those of the 
fully calibrated spectra, for which the dispersion relations 
should be accurate to much better than 1 km s “ 1 (Robinson et 
al. 1992). However, other sources of uncertainty, such as asym- 
metries in the aberrated point-spread function (Gilliland & 
Hulbert 1993) and possible inaccuracies in target centering, 
may produce larger errors than the uncertainties associated 
with the dispersion relations. For our moderate-resolution 
spectra, we estimate the accuracy of our wavelength cali- 
brations to be ±3 and ±5 km s 1 for the fully calibrated and 
SPYBAL-calibrated spectra, respectively. The wavelength cali- 
bration of the echelle spectra should be good to within ± 1.5 
km s ' 1 (Linsky et al. 1 993). 

Most of the spectra were taken in a series of at least two data 
readouts. Breaking up the observations in this fashion permit- 
ted us to check the individual readouts for systematic line shifts 
due to changes in the Earth’s magnetic field and/or instrumen- 
tal thermal drifts. We also could search for flux variations in 
the emission lines, which is especially useful for HR 1099, 
which flares frequently. However, no substantial line shifts or 
line flux variations were evident during any of the observations 
listed in Table 1. The G160M spectra of the 1205-1243 A 
region for Procyon and HR 1099 were obtained using the 
FP-SPLIT mode, in which each observation consists of numer- 
ous readouts with the spectrum shifted on the diode array 
between readouts several times during the observation. This is 
the most accurate way to compensate for sensitivity variations 
among the detector diodes. 

All of the observations in Table 1 are quarter-stepped 
spectra, meaning that each readout of each spectrum consists 
of four separate exposures on the GHRS detector’s array of 
500 science diodes, and the spectrum has been shifted between 
each exposure by one-quarter diode. Thus, the spectra were 
oversampled by a factor of 4, and the final spectra are arrays 



TABLE 1 

Summary of G HRS Observations 


Spectra] Spectral Exposure Start 

Range Resolution Orbital Time Time Important Spectral 

Grating Aperture (A) (km s _l ) Date Phase (s) (UT) Features 

Procyon 

G200M LSA 1880-1920 20 1992 Dec 21 ... 135 14:58 Si ill], C in] 

G160M LSA 1623-1658 20 1992 Dec 21 ... 592 15:56 Hen 

G160M LSA 1532-1568 20 1992 Dec 21 ... 377 16:09 Civ 

G160M LSA 1384-1420 20 1992 Dec 21 ... 431 16:20 Si iv, O iv] 

G160M SSA 1205-1243 15 1992 Dec 21 ... 3231 17:31 Si in, H l, O v], N v 

EB-20 SSA 2792-2807 3.5 1992 Dec 21 ... 646 19:27 Mg n 

EB-22 SSA 2593-2605 3.3 1992 Dec 21 ... 646 19:47 Feu 

HR 1099 

GI60M SSA 1205-1243 15 1993 Sep 14 0.24 3916 18:51 Si ill, H l, O v], N v 

G200M LSA 1880-1920 20 1993 Sep 18 1.64 489 19:19 Si in], C in] 

G160M LSA 1637-1672 20 1993 Sep 18 1.65 707 19:32 He II, C I, O ill] 

G160M LSA 1532-1568 20 1993 Sep 18 1.66 816 20:36 Civ 

G160M LSA 1384-1420 20 1993 Sep 18 1.67 1197 20:53 Si iv, O iv] 

EB-22 SSA 2593-2605 3.3 1993 Sep 19 1.73 979 1:18 Fen 

EB-20 SSA 2792-2807 3.5 1993 Sep 21 2.73 1414 21:05 Mg n 
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Fig. 1. — Moderate-resolution spectra of Procyon with line identifications. The 1 384-1420 A spectrum has been smoothed slightly to reveal the weak O iv] lines. 
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Fig. 2. — As Fig. 1, but for HR 1099 


with 2000 elements. (The final FP-SPLIT spectra actually have 
more than 2000 elements because of the shifting of the spec- 
trum on the diode array between readouts.) 

The Procyon and HR 1099 data were processed with the 
1993 September and 1993 December versions of the CALHRS 
calibration software, respectively (Duncan 1992). The sensitiv- 
ities of the GHRS gratings are probably not known to better 
than ± 10%, so we estimate this to be the uncertainty in our 
flux calibrations. The moderate-resolution spectra are dis- 
played in Figures 1 and 2. All obvious emission lines have been 
identified. The strong lines were previously identified and mea- 
sured in WE spectra (Brown & Jordan 1981 ; Simon & Linsky 
1980). We used the solar spectrum (Sandlin et al. 1986) to help 
identify some of the weak lines. The echelle observations of the 
Mg ii h and k lines of Procyon and HR 1099 are displayed in 
Figure 3. For Procyon, the interstellar absorption features 
found in Lya, Mg n, and Fe n have already been analyzed by 
Linsky et al. (1995a). An analysis of similar features seen in the 
spectra of HR 1099 will be presented in a future paper. 

3. ANALYSIS OF PROCYON DATA 

3. 1 . Line Profile Parameters and the Search for 
Broad Components 

The simplest way to characterize a symmetric line profile 
quantitatively is to measure the velocity, flux, and FWHM of 


the line, and one way to estimate these quantities is to fit a 
Gaussian to the line. We have fit single Gaussians to all of 
Procyon’s emission lines using yf minimization to determine 
the best fits (Bevington & Robinson 1992). The results of this 
procedure are summarized in Table 2. In the fitting process, 
each proposed fit was convolved with the instrumental profile 
extracted from Gilliland & Hulbert (1993) before testing the fit, 
in order to correct for instrumental broadening. The rest wave- 
lengths assumed in this paper are the same as those used by 
Linsky et al. (1995b). The velocities listed in Table 2 assume a 
radial velocity for Procyon of —3.5 km s \ which was com- 
puted from the orbital parameters given by Irwin et al. (1992). 
The 1 a random errors in the tabulated parameters have been 
computed following the prescription of Lenz & Ayres (1992), 
and the quality of the Gaussian fits is measured by the reduced 
X 2 statistic listed in the last column of the table. A Gaussian fit 
was not performed for the Lya line because of its non-Gaussian 
shape and strong interstellar absorption. We have estimated a 
flux for this line, however, using the correction for interstellar 
absorption adopted by Linsky et al. (1995a), which suggests 
that 47% of the stellar emission is absorbed by interstellar 
hydrogen. Inspection of the line wings has also provided an 
estimate of the line velocity (see Table 2). A similar procedure 
was used to estimate the velocities and fluxes of the Mg ii h and 
k lines (see § 3.5). The fact that the Mg n and Lya lines are well 
centered in the stellar rest frame suggests that our wavelength 
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Wavelength 

Fig. 3 — Echelle spectra of the Mg n h and k lines of Procyon (to/?) and HR 1099 (bottom). The dotted lines mark the rest wavelengths of the Mg n lines in the 
stellar rest frames. The line locations are marked for both stars of the HR 1099 system. Because of the orbital velocities of the two stars, the G5 star's lines are 
redshifted relative to the K1 star's lines. Besides the Mg ti h and k lines at rest wavelengths of 2802.705 and 2795.528 A, the location of the Mg n subordinate line at 
2797.998 A is indicated. 


calibrations are very accurate. 

Inspection of the profiles of the O v], N v, Si iv, and C iv 
lines revealed a blue-wing excess (see § 3.2), so the blue wing 
was not considered in the single-Gaussian fits to these lines. 
Figure 4 shows that when the blue wings are ignored, these 
lines are well fit by single Gaussians, with the exception of N v 
A1239. According to the yf values in Table 2, all of the single- 
Gaussian fits to Procyon’s emission lines are acceptable within 
the 2 <y criterion except for the fits to the Si in /1207 and N v 
A 1 239 lines. The non-Gaussian nature of the Si m line is prob- 
ably due to appreciable opacity (see § 3.4). The non-Gaussian 
nature of the N v line is different, however, as it is the only line 
that seems to have excess flux in both wings. Opacity broaden- 
ing cannot account for these broad wings for reasonably low 
optical depths, and N v should in any case be less optically 
thick than the Si m, Si iv, and C iv lines. We therefore interpret 
the broad wing emission as being indicative of a turbulent 
component to the N v line, and we fitted the N v line with two 
Gaussians. The results are illustrated in Figure 4 b and Table 2. 
Monte Carlo simulations were used to estimate errors for the 
two-Gaussian fit parameters of N v. Because of the lack of 
excess emission in both wings of the Si iv and C iv lines and 
because these line profiles appear to be Gaussian everywhere 
except in the blue wings, we conclude that these lines do not 
have detectable “broad components” like those found in the 
lines of AU Mic and Capella. This is consistent with the associ- 
ation of the broad components with microflaring since 
Procyon is much less active than AU Mic and Capella. 

The only emission line that seems to have a broad com- 
ponent is the N v /l 1 239 line, which is the highest temperature 


resonance line that we have observed, with a line formation 
temperature of log T = 5.25. This discovery is consistent with 
that of Linsky et al. (1995b), who found that the broad com- 
ponents of Capella’s emission lines are strongest in the higher 
temperature lines. The O v] A 1218 line, which is formed at a 
temperature of about log T = 5.4, does not appear to have a 
broad component, but this is a low optical depth intersystem 
line rather than a resonance line like N v. If the microflares 
that are presumably responsible for the N v line’s broad wing 
emission have high densities, collisional de-excitation could 
conceivably cause the flux emitted in the O v] line to be much 
less than that observed in N v. Since the O v] line is super- 
posed upon the red wing of Lya, we cannot assume a flat 
background as we have for the other lines but must instead use 
a polynomial fit to estimate the “continuum” for the line (see 
Fig. 4a). Because of the uneven continuum, this background 
estimation is more uncertain than for the other lines, meaning 
it is conceivable that the O v] line does in fact have a weak, 
undetected broad component. Linsky et al. (1995b) were 
unable to find evidence for broad wings in Capella’s inter- 
system lines, but various difficulties prevented them from con- 
cluding that their data were inconsistent with the existence of 
broad components in those lines. 

3.2. Line Shifts and Blue-Wing Excesses 
Figure 5 illustrates the blue-wing excesses of the O v], N v. 
Si iv, and C iv lines by displaying the smoothed residuals of the 
fits shown in Figure 4. The velocity of the excess emission is 
about -90 km s -1 , and the emission is most apparent in the 
C iv lines and the stronger Si iv line. The Si in a 1207 line does 
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TABLE 2 

Line Profile Parameters for Procyon 


Ion 

log T* 

(A) 

^men 

(A) 

t> b 

(km s" 1 ) 

/ 

<10 ,3 ) 

FWHM 
(km s~ *) 

xl 

Single-Gaussian Fits 

Si ih c 

4.7 

1206.510 

1206.512 

4 -F 1 

57.7 + 0.5 

81+2 

1.492 

H i d 

4.2 

1215.670 

1215.658 

0± l 

1300 ± 100 



O v] 

5.4 

1218.344 

1218.373 

11 ± 1 

8.4 + 0.2 

81+2 

0.882 

Nv 

5.25 

1238.821 

1238.845 

9 ± l 

10.7 + 0.3 

80 ± 2 

1.557 

Si jv 

4.8 

1393.755 

1393.768 

6 + 1 

20.2 + 0.5 

75 + 2 

1.320 

O iv] 

5.2 

1401.156 

1401.184 

9 ± 2 

2.7 ± 0.2 

49 ± 4 

0.802 

Si iv 

4.8 

1402.770 

1402.779 

5 4- 1 

11.3 ±0.4 

11 ±2 

0.915 

O iv] ( + S iv]) 

5.2 

1404.806 

1404.804 

3 ± 6 

1.3 + 0.2 

81 ± 11 

0.581 

O iv] 

5.2 

1407.386 

1407.410 

9 ± 7 

0.6 + 0.2 

51 ± 13 

0.709 

Ni n 

4.3 

1411.071 

1411.080 

5 ± 3 

1.5 ±0.2 

49 ±4 

0.874 

Fe ii + Ni ii + S l 

4.3 


1412.853 


1.2 + 0.2 

45 + 6 

0.973 

Si ii 

4.3 

1533.432 

1533.416 

0 ± 2 

3.3 ± 0.4 

45 ±4 

0.883 

C iv 

5.0 

1548.202 

1548.232 

9+ 1 

52.8+ 1.0 

79 + 2 

0.694 

C IV 

5.0 

1550.774 

1550.795 

8 ± 1 

28.7 ± 0.7 

77 ± 2 

1.284 

Fe ii 

4.3 

1559.084 

1 559.089 

4±2 

4.9 + 0.6 

37 ±4 

1.392 

Fe ii 

4.3 

1563.788 

1563.781 

2±2 

3.2 ± 0.4 

31 ±4 

1.064 

He ii 

4.0' 

1640.428* 

1640.429 

4+ 1 

7.0 + 0.5 

38 ±2 

1.261 

Mg u d 

4.2 

2795.528 

2795.500 

0 + 1 

2810 + 50 



Mg ii" 

4.2 

2802.705 

2802.677 

0+ 1 

2310 + 50 



Two-Gaussian Fits 

N v 

5.25 

1238.821 

1238.842 

9+ 1 

6.6 + 0.9 

62 + 4 

0.979 


5.25 

1238.821 

1238.836 

7 ± 4 

4.8 ± 0.9 

141 ± 16 

0.979 


* Approximate temperature ofline formation, based on the peak of the line contribution function. 
b Line velocity relative to the star, assuming a stellar radial velocity of —3.5 kms' 1 . 

c Optical depth at line center was a free parameter of this fit (see § 3.4). The FWHM given in the table is the line 
width that would have been observed if the line was not opacity- broadened, so this is the line width that should be 
characteristic of the thermal and turbulent gas motions. 

d The fluxes and velocities of the Lya and Mg n lines were measured directly rather than by a Gaussian fit. 
e The line formation temperature and rest wavelength listed are for the radiative recombination model (see § 3.2). 


not seem to have any excess emission at —90 km s~ \ so 
log T = 4.7 may represent a lower temperature limit for the 
blue-wing excess phenomenon. With the strong, uneven con- 
tinuum under most of Procyon’s emission lines, it is tempting 
to blame the blue-wing excesses on continuum variations 
and/or line blends, but this does not explain why there are no 
red - wing excesses as well. Therefore, we believe the excess emis- 
sion in the blue wings of the transition-region lines is due to 
emission in those lines. 

Despite the blue-wing excesses, the line velocities listed in 
Table 2 demonstrate that the transition-region lines as a whole 
are redshifted with respect to the photosphere. Ayres (1991) 
also found that Procyon’s transition-region lines are red- 
shifted, although his velocities measured from lower quality 
WE data do not agree particularly well with ours. The line 
velocities are plotted versus line formation temperature in 
Figure 6. The line formation temperatures, which are listed in 
Table 2, are based on the peaks of the line contribution func- 
tions (see, e.g., Brown et al. 1984; Linsky et al. 1995b). 
However, the lines with temperatures listed at and below 
log T = 4.7 are actually formed over a fairly broad range of 
temperatures. The line formation temperatures we have 
assumed for these lines may be too high because the density 
and total mass of emitting material should be higher at lower 
temperatures, thereby allowing more emission to be produced 
at lower temperatures than the contribution functions suggest. 
Nevertheless, Figure 6 demonstrates an obvious increase of 
line redshift with temperature. 

Line redshifts have been observed in the UV spectra of many 
stars besides Procyon, including the Sun (Ayres et al. 1983; 


Ayres, Jensen, & Engvold 1988). For the Sun, the redshifts are 
roughly 10 km s -1 near log T — 5.1 and decrease to both 
hotter and cooler temperatures (Doschek, Feldman, & Bohlin 
1976; Athay 1988; Achour et al. 1995). The redshifts of 
Procyon’s lines are similar to those observed on the Sun, 
although we see no evidence for a decrease in line velocity 
above log T = 5.0. It is unlikely that these redshifts represent a 
net mass downflow. Instead, it is widely believed that we are 
seeing the effects of a circulation system in which the down- 
ward leg of the flow pattern produces most, if not all, of the 
transition-region line emission. This could occur if the 
upflowing material is cool but is heated to transition-region 
temperatures before falling back to the star. Observational evi- 
dence for this scenario centers on solar spicules, which are 
believed to consist of cool material (log T ~ 4.0) that has been 
transported from the chromosphere into the corona. Cheng 
(1992) demonstrated quantitatively that it is possible that the 
spicular material is heated to transition-region temperatures 
before descending back toward the chromosphere, thereby cre- 
ating the redshifts observed in transition-region lines. UV 
images of the solar limb in fact show spicular structure at 
transition-region temperatures as well as chromospheric tem- 
peratures (Feldman, Doschek, & Mariska 1979; Cook & 
Brueckner 1991). 

The outwardly directed solar chromospheric jets detected in 
UV lines of Si n and C i by Dere, Bartoe, & Brueckner (1983) 
are probably spectral signatures of the spicules. The absence 
of spatial resolution prevents us from directly determining 
whether chromospheric jets exist on Procyon. Rather than 
being blueshifted, Procyon’s disk-averaged chromospheric 
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Wavelength 

Fig. 4.— Gaussian fits to six of Procyon’s strong transition-region lines. The data are displayed in histogram form, and the dots show the 1 a errors. The fit to the 
N v A 1239 line is a two-Gaussian fit. The parameters for these fits are listed in Table 2. 


lines are either slightly redshifted or show no Doppler shift 
whatsoever. This by no means precludes the existence of spic- 
ular structures on Procyon, however, since the solar chromo- 
spheric jets would not be easily detectable in disk-integrated 
spectra. It is not known whether the mass flows implied by the 
Mg ii h and k line asymmetries seen in the Sun, Procyon (see 
§ 3.5), and other stars are associated with the chromospheric 
jets and/or spicules (Rutten & Uitenbroek 1991). 

Downflowing spicular material is just one possible explana- 
tion for the redshifted transition-region lines. In addition to the 
numerous hypotheses summarized by Brekke (1993), Reale, 
Peres, & Serio (1995) have recently added another, which 
attributes the redshifts to downflowing material produced 
during the evolution of thermal perturbations in the transition 
region. We also call attention to Mullan & Cheng’s (1994) 
acoustically heated model atmosphere for Procyon, in which 
they predict that the average particle velocity in the transition 
region is 2-9 km s' 1 directed toward the photosphere, which 
may be consistent with our observed redshifts for the 
transition-region lines. Mullan & Cheng explain that in their 
model the material lifted upward by acoustic shocks must 
eventually come down, but the upwardly directed shocked 
material should be denser than the downflowing material. The 
predicted net downdrafts can then be seen as a consequence of 


mass conservation, which requires n u v u = n d v d , where n u and v u 
are the time-averaged density and velocity of the upflowing 
material and rt d and v d are the time-averaged density and veloc- 
ity of the downflowing material, respectively. Since n u > n d , v d 
should be greater than v u . However, since the flux in a col- 
lisionally excited line is proportional to n 2 and the upflowing 
material has the higher density, the net downdraft predicted by 
the model does not necessarily produce a redshifted emission 
line. 

McClymont & Craig (1987), among others, have shown that 
line redshifts can be produced by steady flows through a mag- 
netic loop that is heated predominantly in one leg. In their 
model, the portion of the loop that contains downflowing 
material at transition-region temperatures is much longer than 
the corresponding section that contains upflowing material, 
which results in a much larger emission measure for the down- 
flowing leg of the loop. Perhaps the blue-wing excesses we 
observe in Procyon’s transition-region lines could be a signa- 
ture of the upflowing leg of such a circulation system. Taken at 
face value, the blue-wing excesses seem to indicate the presence 
of outbound mass flows with velocities of ^90 km s' 1 . Such 
flows would be supersonic since the sound speed at transition- 
region temperatures is 30-60 kms" 1 . Outflows such as this are 
occasionally seen on the Sun during transition-region explo- 







768 


WOOD ET AL. 


Vol. 458 



-150 


100 -50 0 50 

Velocity (krn s' 1 ) 


100 


Fig. 5. — Smoothed residuals of the fits in Fig. 4. The horizontal dotted lines 
are the zero flux levels for the six lines. The residuals at the line centers < - 45 
km s~ 1 < t; < 45 km s _ l ) have been set to zero to focus attention on the wings. 
To permit a more equitable comparison of the six emission lines, they were 
normalized to have the same integrated flux before the residuals were com- 
puted. ProcyorTs transition-region lines all seem to have excess flux at about 
— 90 km s 1 ( dashed line). There does not seem to be a corresponding flux 
excess in the red wing. 


sive events (see, e.g., the “jets” discussed by Brueckner & 
Bartoe 1983). However, the blue-wing excesses of Procyon’s 
disk- and time-averaged transition-region line profiles seem to 
be localized to the blue wings of the emission lines, whereas a 
turbulent phenomenon like the explosive events should affect 
the entire profile and presumably produce emission in the red 
wing as well, as discussed in § 3.1. Therefore, we are not sure if 
there is a true solar analog for Procyon’s blue-wing excess 
phenomenon. 

Two very different line formation mechanisms have been 
proposed for the He n 21640 multiplet, which imply very differ- 
ent line formation temperatures (see, e.g., Feldman et al. 1975). 
The collisional excitation model predicts that He n is formed at 
log T ~ 5.0, whereas the radiative recombination model pre- 
dicts line formation at lower temperatures near log T ~ 4.0. 


12 r 


jx 6 


Mg n M 


1 0 A . 5 5.0 5.5 

I-«R 'I' 

Fig. 6. — Line redshifts for Procyon as a function of line formation 
temperature. 


For Procyon, the width of the He ii line is similar to that of the 
Fe ii, Si ii, and Ni n lines, and if we use the centroid wavelength 
of the radiative recombination model (1640.428 A) as the rest 
wavelength of the line, the resulting line velocity, v = +4 km 
s“ l , is also consistent with those of the Fe ii. Si ii, and Ni ii 
lines. If we had instead used the centroid wavelength of the 
collisional excitation model (1640.438 A), the line velocity 
would have been 1 km s -1 , which is not consistent with the 
velocities of the lines formed at log T ~ 5.0. Therefore, we 
believe that the He ii multiplet is mostly formed at chromo- 
spheric temperatures on Procyon. Wahlstrom & Carlsson 
(1994) demonstrated that this is also the case on the Sun. In the 
radiative recombination model, energetic radiation from the 
corona ionizes helium atoms, which then produce the He n 
21640 emission through recombination. In the bandpass of the 
Einstein IPC, the X-ray flux of Procyon is f lPC = 9 x 10“ 12 
ergs cm~ 2 s -1 (Wood et al. 1994), which means that / ipc //hcii = 
13 for Procyon. Hartmann, Dupree, & Raymond (1982) pre- 
dicted that ApcZ/hci ~ 30 for the radiative recombination 
model of He ii line formation. This prediction is too high by 
about a factor of 2, but the uncertainty in the theoretical 
/ipc//hcii ™tio is surely quite large. 

3.3. Electron Density of the Transition Region 

The electron density is one of the fundamental physical 
parameters that must be measured or estimated before a model 
of a stellar atmosphere can be calculated. Spectral density 
diagnostics are therefore of crucial importance. In UV spectra, 
all such diagnostics rely on density-sensitive line ratios, which 
include optically thin intersystem lines such as the O v] 21218 
and O iv] 21400 multiplet lines we have observed in our 
Procyon spectra. Unfortunately, the continuum in the 1880- 
1920 A region proved to be too strong to detect the Si hi] 
2 1 892 and C hi] 2 1 909 lines (see Fig. 1 ). 

The 1395-1410 A spectral region contains four O iv] lines 
that were detected in the GHRS spectra of Capella. The flux 
ratios of these O iv] lines are the most accurate diagnostics 
available in our spectra for determining transition-region den- 
sities (Linsky et al. 1995b). The strongest O iv] line at 1401.156 
A is obviously detected in our spectra, but the weaker lines are 
difficult to see. We have smoothed the 1380-1420 A spectrum 
in Figure 1 to try to reveal these lines. There appears to be a 
feature near 1400 A that could be the O iv] 21399.8 line, but 
the large width of the feature suggests that it is likely blended, 
so we have not listed this feature in Table 2. For similar 
reasons the feature at 1406 A, which could be the S iv] 21406.1 
line, has also not been listed in Table 2. The possible detection 
of S iv] 21406 means that the weaker S iv] line at 1404.8 A may 
be blended with the O iv] line we have identified at that wave- 
length. However, the calculations of Cook et al. (1995) demon- 
strate that the O iv] line should still account for at least 75% of 
the flux in the feature, even if the feature at 1406 A is entirely 
S iv] 21406, so we have identified the line as O iv] 21405 and 
listed the flux of this line in Table 2. The O iv] line at 1407.4 A 
is visible, but it is barely a 3 o detection. Note that the O ivl 
21399.8 line should have the same flux as the line at 1407.4 A 
(Cook et al. 1995). 

The electron density can now be estimated from the 
observed flux ratios of the O iv] lines listed in Table 2 and an 
atomic model of O iv] (Cook et al. 1995). Unfortunately, the 
errors in the fluxes of the weaker O iv] lines are high because of 
the possible S iv] blend in the 1405 A line and low signal- 
to-noise ratios in all the lines. The flux ratios seem to indicate 
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densities in the neighborhood of 10 10 cm" 3 . However, this 
density is not too far from the lower limit where the flux ratio 
versus density relations flatten, so when the large errors in the 
flux ratios are taken into consideration we find we can only 
quote an upper limit for the electron density of log n e < 10.3. 
This provides an upper limit for the transition-region pressure 
(in units of cm - 3 K) of log P e < 15.5. 

Using Fe xix lines observed by EUVE , Schmitt, Haisch, & 
Drake (1994) measured a coronal density for Procyon of 
log n e = 9.70 ±0.15 at a temperature of log T = 6.3. The 
implied coronal pressure is log P e = 16.00 + 0.15, which is sig- 
nificantly higher than our upper limit for the pressure in 
Procyon’s transition region. The interpretation of WE and 
X-ray data presented by Jordan et al. (1986) suggested a similar 
pressure difference between the transition region and corona. 
This pressure difference is inconsistent with a simple, spher- 
ically symmetric, one-dimensional atmosphere in hydrostatic 
equilibrium and instead suggests that Procyon’s transition- 
region emission is not uniformly distributed on the stellar 
surface. Jordan et al. (1986) proposed that this could imply 
supergranulation structures on Procyon similar to those seen 
on the Sun, where most of the emission is produced at the cell 
boundaries. Feldman & Laming (1994) showed that, in general, 
no physical connection appears to exist between the solar tran- 
sition region and the solar corona. If this is true on Procyon as 
well, then the pressures in the transition region and corona 
should not be expected to be consistent with hydrostatic 
equilibrium. 


vective motions should be more vigorous and should penetrate 
to higher altitudes on Procyon than on the solar-like G2 V star 
a Cen A. These properties of convection on Procyon might 
explain why Procyon has broader emission lines than the Sun. 

These estimates of the nonthermal velocities are accurate 
only if the lines are not significantly opacity-broadened. Only 
the Lya and Si hi A 1207 lines of Procyon have profiles that are 
noticeably affected by opacity (see below). The stronger C iv 
line should have twice the opacity of the weaker line, so if the 
lines were optically thick the stronger line should be broader. 
The excellent agreement between the widths of the stronger 
and weaker C iv lines suggests that either these lines are not 
optically thick or that macroturbulence is masking the opacity 
effects. The same is true for the two Si iv lines, the widths of 
which are also in agreement. On the other hand, the C iv 
Al 548/2 1551 and Si iv A1394/A1403 flux ratios suggest the pres- 
ence of some detectable opacity in these lines. Since the ratio of 
collisional excitation rates for the C iv and Si iv multiplets is 
2.0 and these lines do not share a common upper level, the flux 
ratios should be 2.0 if the C iv and Si iv lines are optically thin 
or effectively thin. Instead, we find that these ratios are close to 
1.8, similar to the ratios found for AU Mic (Linsky & Wood 
1994) and Capella (Linsky et al. 1995b) but different from the 
Sun’s line ratios, which are close to 2.0 (Dere et al. 1 987). 

As previously stated, the Si in A 1207 line has a profile that we 
believe is opacity-broadened. Figure la shows a single- 
Gaussian fit to the Si in line. The Gaussian is not a particularly 
good fit (xl = 2.201), because the top of the observed profile is 


3.4. Line Widths and Opacities 

The weighted average width of the emission lines formed at 
temperatures above log T — 4.6 is FWHM = 77.6 + 2.0 km 
s" 1 , where the quoted error is the standard deviation of the 
mean. The O iv] A 1401 and A 1407 lines are the only transition- 
region lines that have widths that seem discrepant, as these 
lines are significantly narrower than the other transition- 
region lines. It is unclear exactly why the O iv] line widths are 
discrepant, but we note that the width of the O iv] lines is 
similar to that of the narrow component of the N v line, which 
is formed at essentially the same temperature as O iv]. 
Perhaps O iv] emission in the broad component is inhibited 
by high densities, as is suggested in § 3.1 for the O v] line. If 
the O iv] lines are ignored, the mean line width becomes 
FWHM = 78.7 ± 0.9 km s" 1 . The chromospheric Ni ii, Fe n. 
Si ii, and He ii lines are much narrower, with an average width 
of FWHM = 39.6 ± 2.3 km s" 1 . Rotational broadening makes 
little contribution to these widths since v sin i is only 3.6 km 
s" 1 (Gray 1982). The nonthermal velocities present in the tran- 
sition region and chromosphere can be estimated for optically 
thin lines using the equation (in cgs units) 



= 3.08 x 10“ 21 




( 1 ) 


where AA is the FWHM, T is the line formation temperature, 
m, is the ion mass, and c is the “most probable” nonthermal 
speed. Using this equation, we find that ^ = 46.1 ± 0.5 km 
s" 1 for the transition-region lines (excluding O iv]) and 
£ = 22.7 ± 1.6 km s“ l for the chromospheric lines. These 
velocities are significantly higher than the solar values: C ~ 25 
km s“ 1 for solar transition-region lines and £ % 15 km s ' 1 for 
solar chromospheric lines (Dere & Mason 1993). The hydrody- 
namic models of Nordlund & Dravins (1990) suggest that con- 
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Fig. 7 — (a) Single-Gaussian fit to the Si m /1207 line of Procyon with the 
residuals of the fit below. The data are displayed in histogram form, and the 
dots show the 1 a errors. The xt value is a measure of the quality of the fit. (b) 
Same as (a), but the optical depth has been included as a free parameter of the 
fit (see § 3.4). The yl value is lower than in (a), indicating a much better fit. This 
illustrates that the non-Gaussian nature of the line can be explained by opacity 
broadening. 
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more rounded and less pointed than the Gaussian. To demon- 
strate that the non-Gaussian shape of the Si iii line can be 
explained by opacity broadening, we have used the following 
technique to perform another fit to the line in which the line- 
center optical depth is a fourth free parameter. 

The formal solution to the radiative transfer equation is 

/ v (t = 0, n) = - | SMe-'^dt . (2) 

P Jo 

If (j) v is the line profile function, the optical depth at line center 
is x 0 = (<£ 0 /$ v )t v . In order to consider the optical depth as a free 
parameter of our profile-fitting routine, we first had to assume 
a functional form for the source function S v (f). We chose to use 
an Eddington-Barbier relation of the form S v (t ) = a + bt with 
a — 1 and b = 4> 0 /<t> v . Equation (2) can be solved numerically 
by using a quadrature sum (Mihalas 1978). For a one-point 
Gaussian quadrature, g = 1/3 1/2 . With these assumptions, 
equation (2) becomes 


/v i'+tvy v * ,dt ~ 

This equation can be solved analytically to yield 

'■ w = ' v [ i + 7 !£ 


(3) 


(4) 


where N is a factor we use to normalize the computed profile to 
the observed profile. This is the equation we used to obtain the 
fit in Figure lb. 

The fit shown in Figure lb is a significant improvement over 
that shown in Figure la. This is not an acceptable 2 o fit, 
however, mostly due to excess emission on the red shoulder of 
the line, which may represent a real asymmetry in the line 
profile. The Si hi line width listed in Table 2, FWHM = 81 km 
s“\ is the width that would have been observed if there were 
no opacity broadening. For comparison, the width of the 
Gaussian in Figure la is FWHM = 92 km s l , illustrating the 
degree to which opacity has broadened the line. Removal of 
opacity broadening leads to an intrinsic width (81 km s~ i ) that 
is consistent with the widths of the other transition-region 
lines. Despite the uncertainties present in our analysis, the fit in 
Figure lb illustrates that the non-Gaussian nature of the Si iii 
line can be explained by opacity broadening. 

The optical depth we derive for Si iii A 1207 is x 0 = 0.7, but 
this optical depth is highly uncertain. The Eddington-Barbier 
relation we used to model the source function is more appro- 
priate for emission in an extended atmosphere. However, a 
slab geometry is more appropriate for transition-region lines, 
since they are emitted in a relatively thin atmospheric layer. 
The source function of a transition-region line should increase 
with depth initially but should reach a maximum near the 
center of the emission layer and then decline until the bottom 
of the layer is reached. We experimented with schematic source 
functions that have this property and found that these source 
functions suggested somewhat higher optical depths for the 
Si in line. 

Another possible problem with equation (4) is that it implic- 
itly assumes that the Doppler broadening is due entirely to 
microturbulence rather than macroturbulence. The line 
opacity is dependent on the former but not on the latter. 


Macroturbulence can be included as a fit parameter by con- 
volving the profile given by equation (4) with a Gaussian, 
where the width of the Gaussian is the fifth parameter of the fit. 
This assumes, however, that the macroturbulence does in fact 
have a Gaussian velocity distribution, which is by no means 
certain. We tried fitting the Si m line with such a profile model 
and found that the inclusion of macroturbulence did not 
improve the quality of the fit. Since macroturbulence gives a 
profile a more Gaussian shape, it can hide the effects that 
opacity has on the line profile. Therefore, if macroturbulence is 
present in the Si iii line, the line’s optical depth will be larger 
than the x Q = 0.7 value we derived. At the very least, the x 0 — 
0.7 result should be a solid lower limit to the actual optical 
depth, since the use of an accurate source function and a cor- 
rection for any macroturbulence present would result in a 
larger derived optical depth. 

The optical depth of a line and the line’s surface flux are 
related by the equation 

t 0 9.7 x 10 ^/A 3 T 3/2 / 1.44 x 10 8 

F = SIMP, eXP [-W- 

where AA is the FWHM of the line (after correcting for opacity 
broadening) and P e is the electron pressure in units of cm ~ 3 K. 
Both A and AA are in angstroms. Brown & Jordan (1981) 
derived this equation in terms of a Gaunt factor, but we 
express it in terms of the statistical weight of the lower atomic 
level g , the oscillator absorption strength /, and the collision 
strength Q. The line fluxes in Table 2 can be converted to 
stellar surface fluxes by using the stellar radius and distance 
listed for Procyon in § 1, although one must make the assump- 
tion that the emission is produced uniformly over the stellar 
surface. We can then solve for the electron pressure in equation 
(5) since we have an estimate of the optical depth of the Si iii 
A 1206 line and all other quantities are known for that line. For 
the Si ill line, Q = 7.5 (Dufton & Kingston 1989), g = 1, and 
/= 1.669 (Morton 1991). We therefore find that P e = 7.0 
x 10 14 cm -3 K at log T — 4.7, where the Si iii line is formed. If 
x 0 = 0.7 is considered to be only a lower limit, then the pres- 
sure is an upper limit (log P e < 14.8) that is significantly lower 
than the upper limit for the pressure derived in § 3.3 by using 
the O iv] lines (log P e < 15.5). 

3.5. The Mg ii h and k Lines 

The Mg ii h and k lines with rest wavelengths in air of 
2802.705 and 2795.528 A, respectively, are shown in Figure 3. 
For Procyon, the chromospheric Mg ii h and k emission fea- 
tures are both located at the center of very broad absorption 
lines produced by photospheric Mg ii absorption. The much 
narrower chromospheric emission lines at the bottom of these 
absorption troughs are self-reversed because of the large 
optical depth of the Mg n-emitting region. Brown & Jordan 
(1981) estimated an optical depth of order x 0 = 5000 for the k 
line, based on IUE observations and the atmospheric model of 
Ayres, Linsky, & Shine (1974). Very narrow interstellar absorp- 
tion features are seen near the center of both Mg ii lines. 
Thanks to the high resolution of the echelle-B grating, these 
interstellar features are easily separable from the central 
reversals of the stellar emission lines, and an analysis of the 
interstellar absorption has been made by Linsky et al. (1995a). 

The blue emission peaks of the self-reversed Mg ii lines are 
significantly stronger than the red peaks. A similar asymmetry 
is seen for disk-averaged solar Mg ii lines (Donnelly, White, & 
Livingston 1994). Solar Ca n H and K lines, which are similar 




No. 2, 1996 


GHRS OBSERVATIONS OF PROCYON AND HR 1099 


771 


to the Mg ii h and k lines in many respects and show the same 
asymmetry, have been studied extensively with ground-based 
observations. Spatially resolved solar observations of the Ca ii 
lines generally do not resemble the disk-averaged profiles. In 
active regions, the emission is strong but relatively symmetric. 
Asymmetric line profiles with stronger blue peaks are more 
typically observed in quiet regions of the chromosphere, where 
there are many small, bright points of Ca ii emission scattered 
over the solar disk. Most, but not all, of the Ca ii profiles 
observed from these “cell grains” show emission almost 
entirely in the blue peak. Rutten & Uitenbroek (1991) gave a 
thorough review of the observations and models of the solar 
Ca ii H and K cell grains. The cell grains typically have life- 
times of a few minutes and often exhibit an oscillatory behav- 
ior that appears to originate in the photosphere. It seems 
likely, therefore, that the stronger Ca n and Mg ii blue peaks 
seen in the disk-averaged spectra of both Procyon and the Sun 
are signatures of hydrodynamical processes in the chromo- 
spheres of these two stars. 

Since the Mg ii line profiles are controlled by opacity and 
atmospheric structure rather than by turbulence, we do not try 


to fit Gaussians to Procyon’s Mg n lines, but we do estimate 
velocities and fluxes for these lines (see Table 2). Inspection of 
the sides of the emission lines reveals no significant Doppler 
shifts. The Mg ii subordinate line at 2798.0 A, which is seen in 
absorption, also appears to be centered in the stellar rest frame. 
The interstellar absorption features were removed before com- 
puting line fluxes. The flux of the k line was measured between 
the k x minima at roughly 2794.8 and 2796.2 A, and the flux of 
the h line was measured between the h { minima near 2802.0 
and 2803.2 A. 

When measuring the Mg ii fluxes, one can assume a flat 
continuum extrapolated between the k { and minima (at 
about 9.5 x 10" 11 ergs cm -2 s _1 A" 1 and 1.15 x 10 -10 ergs 
cm -2 s" 1 A -1 , respectively), or one can simply assume that 
there is no photospheric continuum under the chromospheric 
emission lines. We used the latter assumption in the computa- 
tion of the Mg ii fluxes listed in Table 3 because the Mg ii ions 
in the photosphere and/or lower chromosphere should absorb 
essentially all of the continuum photons near line center, where 
the chromospheric emission lines are located (Linsky & Ayres 
1978). If we had instead extrapolated a flat continuum between 


TABLE 3 


Line Profile Parameters for HR 1099 


Ion 

log 7“ 

^re*t 

(A) 

(A) 

r b 

(km s“ l ) 

/ 

(10~ 13 ) 

FWHM 

(km s 1 ) 

xi 

Single-Gaussian Fits 

Si in 

4.7 

1206.510 

1206.611 

-9 + 2 

8.7 + 0.3 

167 ± 5 

1.105 

H i c 

4.2 

1215.670 

1215.821 

3 + 3 

1200 + 200 



O v] 

5.4 

1218.344 

1218.437 

-11 ± 12 

1.8 ±0.3 

160+ 22 

1.121 

N v 

5.25 

1238.821 

1238.957 

-1 ± 3 

2.8 ±0.1 

124 + 5 

0.978 

Si iv 

4.8 

1393.755 

1393.520 

8 ± 2 

5.5 ±0.1 

179 ± 4 

1.215 

O iv] 

5.2 

1399.774 

1399.507 

1 ± 13 

0.15 + 0.02 

201 + 24 

1.048 

O iv] 

5.2 

1401.156 

1400.889 

1 + 6 

0.76 ± 0.06 

201 ± 24 

1.048 

Si iv 

4.8 

1402.770 

1402.516 

4 ± 3 

2.7 + 0.1 

175 ± 4 

1.048 

O iv] 

5.2 

1404.806 

1404.539 

I + 13 

0.27 + 0.04 

201 ± 24 

1.048 

O iv] 

5.2 

1407.386 

1407.119 

1 + 13 

0.16 + 0.02 

201 + 24 

1.048 

Si ii 

4.3 

1533.432 

1533.208 

13 + 5 

0.70 ± 0.06 

147 ± 10 

0.962 

Civ 

5.0 

1548.202 

1547.949 

8+1 

26.6 + 0.4 

194 + 2 

2.677 

C iv 

5.0 

1550.774 

1550.521 

8± 1 

13.6 ± 0.3 

178 ± 2 

2.677 

C i 

3.8 

1561 mult. 



2.8 + 0.1 



He ii 

4.0 d 

1640.428“ 

1640.155 

5 ± 1 

19.4 ± 0.3 

154 + 2 

2.564 

Ci' 

3.8 

1656.267 

1656.002 

7+ 1 

2.2 + 0.1 

49+7 

1.375 

C i* 

3.8 

1657 mult. 



10.5 + 0.2 



O in] 

4.9 

1666.153 

1665.895 

9 + 8 

1.2 ±0.1 

209 ± 14 

0.956 

Al ii 

4.3 

1670.787 

1670.560 

14 ±4 

2.5 ±0.1 

185 ± 7 

1.156 

Si m] 

4.6 

1892.030 

1891.693 

0 + 4 

2.5 + 0.1 

160 + 6 

1.135 

Si 

3.8 

1900.286 

1899.940 

— 2 ± 5 

0.50 ± 0.09 

54 ± 8 

1.047 

C 111] 

4.8 

1908.734 

1908.231 

-26 ±4 

2.8 + 0.1 

203 ± 8 

1.233 

Si 

3.8 

1914.698 

1914.392 

5 + 6 

0.30 ± 0.08 

53 ± 11 

1.024 

Two-Gaussian Fits 

Si iv 

4.8 

1393.755 

1393.513 

6 + 2 

2.5 + 0.8 

117 + 20 

1.042 


4.8 

1393.755 

1393.513 

6 + 2 

3.1 ±0.8 

240 + 36 

1.042 

Civ 

5.0 

1548.202 

1547.965 

11 + I 

9.6 ± 0.6 

82 ± 5 

1.137 


5.0 

1548.202 

1547.927 

3 ± 2 

18.4 + 0.7 

250 + 4 

1.137 

C iv 

5.0 

1550.774 

1550.537 

11 + 1 

4.8 + 0.3 

82+5 

1.137 


5.0 

1550.774 

1550.499 

3 ± 2 

9.2 ± 0.3 

250 + 4 

1.137 

He ii 

4.0“ 

1640.428 d 

1640.185 

11 + 1 

7.0 ± 0.6 

52 + 7 

1.120 


4.0“ 

1640.428“ 

1640.119 

-2 + 2 

13.4 ±0.5 

204 +5 

1.120 

C ill] H- ?.... 

4.8 

1908.734 

1908.396 

0 

2.1 ± 0.2 

187 ± 15 

1.046 


? 

? 

1907.795 

7 

0.7 ± 0.2 

101 ± 26 

1.046 


* Approximate temperature of line formation, based on the peak of the line contribution function. 
b The line velocity relative to the K1 star. 

c The flux and velocity of Lya were measured directly rather than by a Gaussian fit. 

d The line formation temperature and rest wavelength listed are for the radiative recombination model (see 
§ 4.2). 

* The C l A 1656.3 line is the only resolvable member of the / 1 657 multiplet. 
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the k v and hy minima, the fluxes of the h and k lines would have 
been 8.87 x 10~ u and 1.48 x 10" 10 ergs cm -2 s" 1 , respec- 
tively. However, the h and k lines have different profiles if we 
assume such a continuum since the h line would appear to have 
a deeper central reversal and a larger blue/red peak ratio than 
the k line. On the other hand, the h and k line profiles have 
essentially the same shape if one assumes a zero continuum 
level for both lines. This is empirical evidence in favor of this 
assumption, since the h and k lines of most stars, including 
Capella (Linsky et al. 1995b) and HR 1099, appear to be nearly 
identical if the flux and width of the h line are slightly 
increased. 

4. ANALYSIS OF HR 1099 DATA 

4.1. Line Profile Parameters and the Search for 
Broad Components 

The first step in analyzing the HR 1099 data was to deter- 
mine which of HR 1099’s two stars is responsible for the 
binary’s emission lines. By analyzing IUE spectra at different 
orbital phases, Ayres & Linsky (1982) found that the K1 star 
accounts for most of the emission. Our data not only confirm 
this finding but also suggest that the G5 star’s flux contribution 
is truly negligible in many of the strongest emission lines 
observed with the moderate-resolution gratings (e.g., C iv, 
Si iv, and He ii). Both the near-symmetric appearance of these 
lines and the good agreement of the line velocities with the rest 
frame of the K1 star support this conclusion. The only lines in 
which the presence of the G5 star is obvious are the Mg ii h 
and k lines, but even for those lines the G5 star accounts for 
only ~7.5% of the line flux (as will be discussed in § 4.5). For 
the C iv, Si iv, and He ii lines, we believe the G5 star’s contri- 
bution must be only ~5% at most. Any larger contribution 
should have been detected because our observations were 
taken at or near quadrature when the velocity separation 
between the two stars is near the maximum of 1 1 1.1 km s~ l . 

We fit Gaussians to all of HR 1099’s emission lines, as we did 
for the Procyon data, and Table 3 lists the results. We cor- 
rected for the instrumental broadening and also corrected for 
rotational broadening by convolving proposed Gaussian fits 
with both the instrumental and rotational broadening profiles 
before testing the fits. This required knowledge of v sin i for the 
stars in the HR 1099 system; v sin i = 38 km s _1 for the K1 
star, and v sin i = 13 km s~ 1 for the G5 star (Vogt & Penrod 
1983). (This latter value was required only for the analysis of 
the Mg ii lines in § 4.5.) An example of the instrumental and 
rotational broadening profiles is shown in Figure 8 for the 
G160M spectra. 

As was the case for Procyon, the Lya line was not fit with a 
Gaussian, but its velocity was estimated from parts of the line 
profile uncontaminated by interstellar absorption. In another 
paper, this emission line and the superimposed interstellar H i 
and D i absorption features will be analyzed in the same way as 
the Lya lines of Capella and Procyon (Linsky et al. 1993, 
1995a), and this analysis will provide an estimate for the 
amount of stellar Lya flux that has been removed by the inter- 
stellar absorption. For now, we use a cruder estimate based on 
the Capella and Procyon work, in which it was found that 
interstellar absorption removed about 60% and 50% of the 
Lya line flux for Capella and Procyon, respectively. We expect 
a similar value for HR 1099 and have therefore assumed a 
similar correction to estimate the intrinsic stellar Lya flux in 
Table 3. The only C i line that is resolvable from either of the 



Fig. 8. — Instrumental profile for the G160M grating (Gilliland & Hulbert 
1993). The profile is shown as a smooth dotted line, and the rotational 
broadening profile for HR 1099’s K1 star (assuming v sin i = 38 km s" 1 ) is the 
dotted histogram. The final line-spread function (solid histogram ) for the 
G160M observations of HR 1099’s K.1 star is the convolution of these two 
profiles. 

C i multiplets at 1561 and 1657 A is the C i Al 656.267 line. The 
line parameters for this line are listed in Table 3, and the fluxes 
for the C i A1561 and A1657 multiplets are also listed. 

Based on the xl values in Table 3, the single-Gaussian fits to 
five of HR 1099’s strongest emission lines are not acceptable 
within the 2 o criterion. Four of these lines (Si iv A 1394, C iv 
AA1548, 1551, and He n A1640) have broad wings that we inter- 
pret to be indicative of “ broad components ” similar to those 
found for AU Mic and Capella. These four lines have therefore 
been fit with two Gaussians, and Monte Carlo simulations 
were used to estimate the 1 o errors in the fit parameters. The 
results are shown in Figure 9 and in Table 3. The fifth non- 
Gaussian line, C ill] A1909, was also fit with two Gaussians (see 
Fig. 9 d and Table 3), but the interpretation of the C hi] line 
profile is completely different from that of the four line profiles 
modeled with narrow and broad Gaussian components (see 
§ 4.2). 

The C iv A1548/A1551 flux ratio is exactly 2.0, the optically 
and effectively thin limit, and the two C iv lines have almost 
exactly the same profile. Therefore, we fit the two C iv lines 
simultaneously and forced the A 1548/A 1551 flux ratio to be 2.0 
for both the broad and narrow components, and we forced the 
widths and velocities of the two components to be the same for 
both C iv lines. The single-Gaussian fit to the Si iv A 1394 line 
falls just short of being a 2 a fit. To reduce the number of free 
parameters for the two-Gaussian fit to that line, we required 
the broad and narrow components to have the same velocity. 
When we tried to do that for the C iv and He ii lines (for which 
the single-Gaussian fits are not even close to being acceptable), 
we could not obtain 2 a fits. Therefore, the velocities of the 
narrow and broad components were allowed to be different for 
those lines. 

We cannot rule out the possibility that HR 1099’s weaker 
lines also have broad components that were not detected sta- 
tistically because of low signal-to-noise ratio. Figure 10 shows 
that when these weaker lines are smoothed, many of them do 
in fact seem to show evidence for non-Gaussian profiles. Com- 
pared to the single-Gaussian fits, the smoothed profiles of N v 
A1239, Al II A1671, and Si iii] A1892 seem to have excess flux at 




1638 1639 1640 1641 1642 1906 1907 1908 1909 1910 

Wavelength 

Fig. 9.— Two-Gaussian fits to four lines of HR 1099 that were not well fit with single Gaussians. The data are displayed in histogram form, and the dots show the 
1 a errors. For each fit, the Gaussian components are shown as thin solid lines, and the thick solid line is the convolution of the sum of these Gaussians with both the 
instrumental and rotational broadening profiles (see Fig. 8). The Si iv, C iv, and He n lines were fit with a narrow component and a broad component. The stronger 
component of the fit to the C in] feature was forced to be centered in the rest frame of the star. We believe that this component represents the actual “quiescent” 
Cm] emission from HR 1099. The weaker component may represent blueshifted C m] emission produced by some transient event, or it may represent a blend. 
The parameters for all these fils are listed in Table 3. 
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Wavelength 

Fig. 10. — Smoothed profiles of four of HR 1099’s UV emission lines and the single-Gaussian fits to the lines. The data are displayed in histogram form, and the 
dots show the 1 a errors. Although these fits are statistically acceptable, smoothing the profiles seems to reveal some non-Gaussian behavior in all these lines. N v, 
A1 H, and Si ill] show at least some evidence for excess flux at line center and in the line wings, which may indicate the presence of broad components such as those 
found for Si iv, C iv, and He n (see Fig. 9). The O in] line has a strange, boxlike appearance. 
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line center and in the line wings, which may indicate a shape 
like that of the C iv and He n lines. We do not feel justified in 
fitting these lines with two Gaussians, however, since the 
single-Gaussian fits are statistically acceptable. We do not 
have an explanation for the strange, boxlike appearance of the 
O hi] A 1666 line in Figure 10 b. 

The broad components of the transition-region lines of AU 
Mic and Capella account for roughly 40% of the line flux 
(Linsky & Wood 1994; Linsky et al. 1995b). In contrast, the 
broad components of HR 1099’s lines actually dominate the 
line flux. The weighted average contribution of the broad com- 
ponent to the total line flux is 65.6% ± 1.0% for the four lines 
in which the broad component was detected. This dominance 
is consistent with the hypothesis that microflares are 
responsible for the broad components because HR 1099’s K1 
star is even more active than AU Mic and Capella’s G1 star on 
the basis of the X-ray and UV emission-line surface fluxes of 
these stars. 

We note that magnetic broadening (i.e., Zeeman splitting) 
and collisional broadening (i.e., Stark and Van der Waals 
broadening) are almost certainly not responsible for the large 
transition-region line widths. A magnetic field of order 10 6 G 
(Brueckner et al. 1988) is necessary to produce lines as broad as 
those of HR 1099, which seems extreme. The line profiles show 
no double-peaked structure that would indicate Zeeman split- 
ting, and even though rotational broadening could help hide 
such structure, the profile of a magnetically broadened 
transition-region line should still look different than the line 
profiles that we observe. Brueckner et al. (1988) pointed out 
that a magnetically broadened C iv A 1551 line should be 
broader than the C iv A 1548 line, which is not observed for HR 
1099 or for any other star that we know of. They also point out 
that for the Stark effect to become detectable densities must be 
at least 10 17 cm -3 , which is much higher than the densities of 
10 10 — 10 1 1 cm -3 typically found in stellar transition regions, 
including that of HR 1099 (see § 4.3 and Byrne et al. 1987). 
Van der Waals broadening is certainly negligible because 
transition-region plasmas are almost completely ionized, 
meaning there are essentially no neutrals available to produce 
this type of line broadening. 

The broad components of HR 1099’s C iv and He n lines are 
significantly blueshifted with respect to the narrow com- 
ponents. This is equivalent to saying that the blue wings of 
these lines are noticably stronger than the red wings. Linsky et 
al. (1995b) noted a similar phenomenon in many of the 
transition-region lines of Capella. The G5 star of HR 1099 
cannot account for this effect since any emission from that star 
would be on the red side of the observed C iv and He ii lines for 
the orbital phase of those observations and would therefore 
produce a stronger red wing rather than a stronger blue wing. 

Individual solar transition-region explosive events produce 
a wide variety of line profiles; some with particularly strong 
blue wings, some with particularly strong red wings, and some 
with roughly symmetric wing emission. However, Dere et al. 
(1989) found that line profiles with dominant blue wings are 
produced more often than profiles with dominant red wings. 
Thus, a disk-averaged solar transition-region line profile 
should have a blueshifted “ broad component,” which is what 
is also observed for Capella and HR 1099. 

4.2. Line Velocities 

The velocities listed in Table 3 are all relative to the rest 
frame of the K1 star, which is given by the equation 


v = — 15.0 + 49.4 sin 2n<f), where (f> is the orbital phase of the 
observation (see Table 1). The rest frame of the G5 star is 
v - — 15.0 — 61.7 sin 2n<f) (Fekel 1983). It is difficult to discuss 
the line velocities of HR 1099 in the same detail as we discussed 
those of Procyon because the uncertainties in the velocities of 
many of HR 1099’s lines are larger than was the case for 
Procyon, because of the larger widths and lower signal-to- 
noise ratios of the lines. Also, the small flux contributions of 
the G5 star may be affecting the measured line velocities. Even 
if the G5 star accounts for only a few percent of the fluxes of 
the emission lines, that is enough of a contribution to produce 
noticeable line shifts at quadrature when our observations 
were made. Using numerical experiments, we estimate that the 
difference in observed line velocity between the two quadra- 
tures should be ^ 1 1 km s ' 1 if the G5 star accounts for 5% of 
the line flux, ~7 km s' 1 if it accounts for 3% of the line flux, 
and ~2 km s' 1 if it accounts for just 1%. Most of the lines 
observed at <j) = 0.24 (Si in A 1207, O v] A 121 8, and N v A 1239) 
are slightly blueshifted, whereas almost all of the lines observed 
near the opposite quadrature (<f> = 0.64-0.67) are redshifted 
relative to the K1 star, which suggests that the G5 star’s flux 
contribution is indeed a factor that must be considered. The 
relatively large blueshift of the Si in A 1207 line and the rela- 
tively large redshifts of the Si ii A 1533 and Al n A 1671 lines 
suggest that the G5 star may be a larger flux contributor to 
these lines than it is to the C iv, Si iv, He n, and N v lines. 

Because the G5 star appears to be at least partially 
responsible for the observed line shifts, the redshifts seen for 
the lines observed at <t> = 0.64-0.67 are actually just upper 
limits for the intrinsic redshifts of the K1 star’s emission lines. 
The K1 star’s transition region lines may have small redshifts 
of -5 km s' 1 , but they certainly do not have redshifts as large 
as those observed on the G1 star of Capella (10-25 km s' 1 ; 
Linsky et al. 1995b). 

For Procyon, the position and width of the He n A 1640 
feature were used to determine whether the line is a col- 
lisionally excited line formed at transition-region temperatures 
or a recombination line formed at chromospheric tempera- 
tures. The situation is more complicated for HR 1099 because 
of the uncertainty in the centroid of the K1 star’s emission 
lines, due to the unknown contribution of the G5 star, and 
because some of HR 1099’s chromospheric lines are just as 
broad as the transition-region lines. The line profiles of the 
He ii, C iv, and the Mg ii h and k lines are very similar for 
reasons that we do not fully understand (see § 4.5). However, 
the width of the He ii line is closer to that of Mg ii than it is to 
the broader C iv lines, suggesting that He ii is a chromospheric 
line such as Mg ii. Therefore, we have used the centroid wave- 
length predicted for the He ii multiplet by the radiative recom- 
bination model as the rest wavelength of the He ii line. The 
average of the two Einstein IPC X-ray fluxes listed for HR 1099 
by Majer et al. (1986) is / IPC = 1.2 x 10' 10 ergs cm' 2 s' 1 , 
which means f lPC /f He h = 59. This is larger than the ratio pre- 
dicted for the radiative recombination model by Hartmann et 
al. (1982) by about a factor of 2. Thus, the corona of HR 1099 
easily produces enough energetic radiation to account for the 
observed He n emission, which further supports the radiative 
recombination model. 

The line velocity of the C m] A 1909 line is -26 km s' 1 , 
according to the single-Gaussian fit (see Table 3). It seems 
unlikely that C hi] would have a velocity so discrepant from 
the velocities of HR 1099’s other lines. The velocities of the 
other three lines identified in the spectral region of C m] (the 
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Si in] A 1892 line and the two S i lifies) are very reasonable, so 
an inaccurate wavelength calibration cannot be responsible. At 
the orbital phase of the observation (</> = 0.64) any emission 
from the G5 star would create an apparent redshift rather than 
a blueshift, so the G5 star cannot be responsible for the dis- 
crepant velocity either. The feature at 1908 A appears to be 
somewhat asymmetric, which may imply a blend of some sort. 
Therefore, we fit the feature with two Gaussians — one rep- 
resenting C in] emission forced to be centeed on the rest frame 
of the K1 star (as the Si in] line appears to be) and the other 
representing the emission blended with C hi]. The results of the 
fit are illustrated in Figure 9 d and Table 3. 

We are not certain what emission line the smaller Gaussian 
in Figure 9 d represents. It could be an unidentified blend, or it 
could represent blueshifted C hi] emission produced by some 
transient event on the K1 star. Although the line formation 
temperature of C in] is somewhat higher than that of Si ill], it 
is still difficult to understand why such an event would not 
have also produced Si in] emission. The observation of the 
1880-1920 A spectral region consisted of two readouts of 
roughly 4 minutes each. There are no large, obvious profile 
differences between the two readouts that would indicate that 
flares or other transient processes are at work here, but the 
individual readouts are very noisy, so it is difficult to say for 
sure. 

4.3. Electron Density of the Transition Region 

Our attempt to derive a transition-region density for 
Procyon in § 3.3 using the O iv] lines was not very successful 
because of the weakness of the lines and the presence of blends. 
The situation is better for HR 1099 because there are no 
obvious blends. In particular, there is no sign of the S iv] 
A1406.1 line, which means we do not have to worry about the 
S iv] A 1404.8 line contaminating the O iv] line at that wave- 
length. The weak O iv] lines are noisy, though, so the 1384— 
1420 A spectral region in Figure 2 has been smoothed 
somewhat to reveal these lines. We fitted all four O iv] lines 
simultaneously (and also the Si iv A 1 403 line that is slightly 
blended with the O iv] lines) and forced the O iv] lines to have 
the same width and velocity to further constrain the fit (see 
Table 3). The fluxes of the weaker O iv] lines of HR 1099 have 
been measured accurately enough to obtain a good density 
estimate, which was not the case for Procyon. We find that the 
fluxes of the A1399.8 and A1407.4 lines are essentially the same, 
which is consistent with the atomic model (Cook et al. 1995). 

Three density-dependent flux ratios can be constructed: 
Ri =/i407//i40i = 0.204 ±0.025, R 2 =/ 1407 //i405 « 0.57 ± 0.10, 
and R 3 =/i 405 //i 40 i = 0.355 ± 0.060. Since the fluxes of the 
>11 407.4 and A 1399.8 lines are the same, as they should be for all 
densities, we use the average of the fluxes of the two lines for 
/i 407 - The electron densities implied by these three ratios are 
log n e - 10.04^;!? for R u log n e = 9.83!g;^ for R 2 , and log 
n e = 9.74±o 23 f° r R 3 ■ These three densities are reasonably 
consistent with each other, and they collectively imply a 
density of log n e = 9.85 ±0.15, which indicates a pressure at 
log T = 5.2 of log P e — 15.05 ±0.15. This pressure is nearly 
identical to pressures found for Capella’s transition region and 
the quiet solar transition region (Linsky et al. 1995b; Dere & 
Mason 1993). We were surprised that HR 1099, Capella, and 
the Sun all have similar transition-region densities and pres- 
sures, considering that the emission-line surface fluxes of HR 
1099 are about 4 times as large as those of Capella and more 
than 100 times larger than those of the Sun. Since the very 
large surface fluxes of HR 1099 do not appear to be caused by 


higher densities and pressures, they must be due to a larger 
filling factor and/or a more radially extended emitting region. 

Using the Si iv AA(1394 + 1403)/C ill] A 1909 flux ratio for 
HR 1099, Byrne et al. (1987) derived a pressure at log T = 4.8 
of log P e = 15.3. The Si iv and C ill] fluxes we have measured 
are similar to the fluxes they measured from their IUE data, so 
we would derive a similar pressure from our data. The pressure 
at log T = 4.8 is slightly higher than the pressure at log T = 
5.2, which is what is expected for an atmosphere in hydrostatic 
equilibrium because of the weight of the overlying atmosphere 
(see, e.g., Brown & Jordan 1981). 

4.4. Line Widths and Opacities 

Unlike Procyon, the widths of HR 1099’s emission lines 
cannot be neatly divided into narrow chromospheric features 
and broader transition-region lines. The only lines that can be 
considered narrow are the weak S 1 AA1900, 1915 lines and the 
C 1 A 1656 line, which have widths of ~50 km s~\ These are 
certainly chromospheric lines, but the Al 11 A 1671 and Si 11 
A 1533 lines, which are also considereed chromospheric, are as 
broad as many of the transition-region lines (see Table 3). It is 
possible that this difference in width among the chromospheric 
lines is due to a dramatic increase in turbulence between the 
lower and upper chromosphere, as the S 1 and C 1 lines are 
probably recombination lines formed at cooler temperatures 
(log T ~ 3.8; Judge 1986, 1988) than the collisionally excited 
Al n and Si 11 lines (log T < 4.3; Brown et al. 1984). If this 
interpretation is correct, then the atmospheric structures of 
both Procyon and HR 1099 show large increases in turbulence 
within a relatively narrow temperature range (but possibly a 
substantial pressure range). Procyon appears to have a factor 
of 2 increase in turbulence somewhere between log T = 4.3 
(represented by the Si 11 A 1533 line) and log T = 4.7 (sampled 
by the Si 111 A 1207 line) whereas HR 1099 shows a factor of 
3 increase in turbulence somewhere between log T = 3.8 and 
log T = 4.3. 

If the Al 11 and Si 11 lines have higher opacities than the S 1 
and C 1 lines, opacity broadening might explain the larger line 
widths rather than greater turbulence. Harper et al. (1995) 
demonstrated that the S 1 A 1 900/A 1915 flux ratio is sensitive to 
the line opacity. Based on their Figure 5, the observed S 1 flux 
ratio for HR 1099 implies an optical depth of t ~ 1 for the 
stronger line. Unfortunately, we have no way of estimating the 
optical depths of the C 1 , Si 11 , and Al 11 lines at this time. Future 
observations of optically thin chromospheric intersystem lines 
such as C 1 ] A1994, Si 11 ] A2350, and Al 11 ] A2669 may further 
our understanding of the physical processes taking place in the 
chromosphere of HR 1099. 

After the S 1 and C 1 lines, the narrowest line in HR 1099’s 
moderate-resolution spectra is N v A 1239, which is signifi- 
cantly narrower than other lines formed at similar tem- 
peratures (e.g., C iv, O iv]). Before trying to explain this 
problem, we point out that the spectral regions shown in 
Figure 2 were not observed simultaneously. Because HR 1099 
is a highly variable object, it is conceivable that this variability 
could contribute to the observed differences in line profile 
parameters, such as the one noted for N v. We have no evi- 
dence that this is what is happening, however, especially since 
the Si hi A 1207 line observed at the same time as N v has a 
width similar to those of most of the observed transition- 
region lines. We will assume, therefore, that we can inter- 
compare line profiles observed in different spectral regions, 
although we cannot completely rule out the possibility of 
variability within our data set. 
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The weighted average width of the lines listed in Table 3 
(excluding the S i and C i lines) is 173 + 19 km s 1 , which 
corresponds to a most probable nonthermal velocity of 
103 + 11 km s 1 . This velocity is an accurate estimate of the 
turbulence in the K1 star’s upper chromosphere and transition 
region if the emission lines are not significantly opacity- 
broadened. With the exception of Lya, the line profiles of HR 
1099 show no obvious signs of opacity effects, as long as one 
assumes that the broad wings of C iv, Si iv, and He u are not 
produced by opacity effects but rather are caused by particu- 
larly turbulent regions of the stellar atmosphere, as we 
assumed in § 4.1. Unfortunately, the Si iii 21207 line is too 
noisy to detect effects of opacity broadening as subtle as those 
seen for Procyon’s Si in line. 

Unlike Procyon, the C iv 21548/21551 and Si iv /U394/ 
21403 flux ratios are exactly 2.0, which is the optically and 
effectively thin limit. The profiles of the two C iv lines and the 
two Si iv lines are identical, which also suggests that the lines 
are optically thin. This does not prove that r 0 < 1 for the C iv 
and Si iv lines since macroturbulence can mask the effects of 
opacity on the line profiles, but there are certainly no obvious 
opacity effects present in the profiles of any HR 1099’s 
transition-region lines. Furthermore, the optically thin inter- 
system lines (O v], O iv], O hi], Si m], and C in]) are just as 
broad as the Si iii, N v. Si iv, and C iv resonance lines, which 
demonstrates that enough turbulence is present in the atmo- 
sphere of the K1 star to account for the broad widths of the 
transition-region lines without any opacity broadening. 

Linsky et al. (1995b) found that the widths ofCapella’s inter- 
system lines were closer to the widths of the narrow com- 
ponents detected in the resonance lines than to the widths of 
broad components. They proposed that the intersystem lines 
are formed predominantly in the narrow-component regions of 
the stellar atmosphere, implying that the density of the broad- 
component regions may be high enough to reduce the inter- 
system line flux substantially relative to the resonance-line flux. 
However, for various reasons, they were unable to demonstrate 
conclusively that the intersystem lines do not indeed have 
broad components. For HR 1099, we are also unable to deter- 
mine conclusively whether the intersystem lines have broad 
components. Some evidence for extended wing emission is 
present in the smoothed Si iii] profile in Figure 10J, but this is 
the only smoothed intersystem line profile that appears similar 
to the resonance line profiles. Unlike Capella, the widths of HR 
1099’s intersystem lines are closer to the widths of the broad 
components than to the widths of the narrow components. The 
O iii], O iv], and C iii] lines are in fact the broadest lines in our 
spectra, based on the single-Gaussian fits. This suggests that 
the less turbulent regions of the stellar atmosphere may have 
pressures high enough to inhibit intersystem line emission and 
that the pressure of log P e = 1 5.05 measured using the O iv] 
line ratios applies only to the broad-component regions. In 
particular, we note that the large discrepancy between the 
widths of the N v and O iv] lines alluded to above suggests 
that these lines are not originating entirely from the same 
regions, despite the fact that the temperatures at which the 
lines are formed are nearly identical. 

4.5. The Mg ii h and k Lines 

The Mg ii h and k line profiles of HR 1099 are displayed in 
Figures 3 and 11. In Figure 11, we have removed the inter- 
stellar absorption features by using polynomial fits to estimate 
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Fig. 11. — The Mg h (bottom) and k (top) line profiles of HR 1099 plotted on 
a heliocentric velocity scale. The 1 o errors in the data are shown as dots. 
Interstellar absorption at roughly +20 km s" 1 has been removed from both 
lines by use of a polynomial fit to estimate the intrinsic stellar emission above 
the interstellar absorption. Three-Gaussian fits were performed for both lines. 
The two components near —60 km s - 1 model the emission from the K1 star, 
and the single component near + 45 km s“ 1 represents the emission from the 
G5 star. The convolution of the sum of the Gaussians with the instrumental 
and rotational broadening profiles is shown as a thick line. The parameters for 
these fits are listed in Table 4. 


the intrinsic stellar emission in the spectral regions obscured by 
the absorption. The Mg ii lines are the only lines we have 
observed in which emission from the G5 star is apparent. We 
performed multi-Gaussian fits to these lines to estimate the 
contributions of the K1 and G5 stars. We found that the K1 
star’s Mg ii profiles could be accurately represented as the sum 
of two Gaussians. A third Gaussian was used to represent the 
G5 star’s emission. The results of the fits are shown in Figure 
1 1 and Table 4. Also listed in Table 4 are the line profile 
parameters of the Mg ii subordinate line at 2798 A. 

The G5 star accounts for only —7.5% of the Mg ii line flux, 
according to the multi-Gaussian fits. However, since the radius 
of the K1 star is 3 times that of the G5 star (R Kl — 3.9 R© and 
Rqj = 1.3 R 0 ; Fekel 1983), the Mg ii surface fluxes of the K1 
star are only ~37% higher than those of the G5 star. Ayres et 
al. (1995) found a power-law slope of 1.77 ± 0.12 when stellar 
C iv surface fluxes were plotted versus Mg ii surface fluxes. If 
such a relation applies to the two HR 1099 stars, then the G5 
star should account for ~6% of HR 1099’s C iv flux, based on 
its 7.5% contribution to the Mg ii lines. This 6% contribution 
may just barely be consistent with the nondetection of the G5 
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TABLE 4 


Mg n Line Profile Parameters for HR 1099 


2 rest 

(A) 

Component 

(A) 

u* 

(km s' 1 ) 

/ 

(I0’ 2 ) 

FWHM 
(km s' 1 ) 

x 1 . 

2795.528 

K1 (narrow) 

2794.959 

3 ± 1 

9.8 ± 0.2 

62 ± 1 

1.494 


K1 (broad) 

2794.995 

7 ± 1 

20.0 ± 0.3 

184 ± l 

1.494 


G5 

2795.945 

-1 ± 1 

2.2 ± 0.1 

62 + 2 

1.494 

2802.705 

K1 (narrow) 

2802.129 

2 4- 1 

7.4 + 0.3 

58 + 2 

1.300 


K1 (broad) 

2802.140 

4 + 1 

13.8 ±0.3 

167 ± l 

1.300 


G5 

2803.144 

1 ± 1 

1.9 + 0.1 

59 ± 2 

1.300 

2797.998 

K1 

2797.414 

1 + 1 

0.35 + 0.01 

34 + 3 

1.395 


G5 

2798.426 

0± 3 

0.026 + 0.006 

26+5 

1.296 


* Line velocity relative to the star responsible for the emission. 


star’s C iv emission, but we believe this value is probably a bit 
large (see § 4.1). 

The multi-Gaussian fits suggest that the K1 star’s Mg ii h 
and k lines are slightly redshifted, with the k line showing the 
larger shift. Similar Mg n redshifts have been observed for 
many other late-type stars (Robinson & Carpenter 1995). The 
fits to HR 1099’s lines also suggest that the broad component is 
more redshifted than the narrow component. It is difficult to be 
certain that these properties of the fits are real, however, 
because of the unquantifiable uncertainties introduced into the 
fits by the assumption that the G5 star’s emission is exactly 
Gaussian and the K1 star’s emission is exactly the sum of two 
Gaussians. The multi-Gaussian fits suggest no Doppler shifts 
for the G5 star’s Mg n lines. 

The two-Gaussian fits to the K1 star’s Mg n lines look strik- 
ingly similar to the fits to the K1 star’s C iv lines (see Fig. 9b). 
This similarity is shown explicitly in Figure 12. Although we 
found it necessary to broaden the Mg n profile by 30% to 
match the width of the C iv line, the agreement between the 
C iv and Mg n profiles in Figure 12 is still remarkable and 



Velocity (km s" 1 ) 


Fig. 1 2. — The C iv (solid line ) and Mg n k ( dotted line) line profiles of the K 1 
star of HR 1099, based on the two-Gaussian fits to these lines. The Mg u line 
has been broadened by 30% to allow a better comparison of the C iv and Mg ii 
line profiles. If it were not for different Doppler shifts in the cores and wings of 
the lines, the two profiles would be essentially identical. The similarity between 
these two profiles is unexpected since the Mg ii profile is not generally sup- 
posed to be regulated by the same turbulent processes that presumably regu- 
iate the C iv profile. 


completely unexpected. The differences that exist between the 
profiles are mostly due to the different Doppler shifts present in 
the cores and wings of the lines— the overall profile shapes are 
nearly identical. As discussed in § 3.5, the Mg ii h and k lines 
are very optically thick lines whose profiles are generally inter- 
preted in terms of opacity effects and atmospheric structure 
whereas the C iv line profiles are thought to be optically thin 
and regulated by turbulent processes. Therefore, the profiles of 
the Mg ii and C iv lines should in general be different, but that 
does not appear to be the case for HR 1099. 

Is the similarity between the Mg n and C iv lines just a 
coincidence, or does it imply that the same physical processes 
are controlling the profiles of both lines? If not a coincidence, 
are the C iv lines behaving like optically thick chromospheric 
lines or are the Mg ii lines behaving like transition-region 
lines? The former case would suggest that the broad wings 
observed for both C iv and Mg ii are produced by opacity 
effects, whereas the latter case would suggest that the broad 
wings are due to particularly turbulent regions of the stellar 
atmosphere. It is important to determine which of these cases is 
correct, because if the broad wings are not caused by turbu- 
lence, then our attempts to explain the broad wings in terms of 
microflaring and to quantify this microflaring by using two- 
Gaussian fits are misguided. 

We have computed theoretical Mg n profiles based on a 
simple, homogeneous, hydrostatic, plane-parallel model atmo- 
sphere of the K1 star to see if radiative transfer calculations 
predict broad line wings similar to those we have observed. 
Following Simon & Linsky (1980), we have based our atmo- 
spheric model on the model for e Eri (K2 V) computed by 
Simon, Kelch, & Linsky (1980). Since the surface gravity of HR 
1099’s K1 star (log g = 3.40; Fekel 1983) is roughly 12.6 times 
lower than that of € Eri, the mass column densities given for 
the € Eri model were increased by a factor of 12.6 to maintain 
the same total pressure as the e Eri model. This scaling 
increases the computed surface fluxes, which is desirable since 
HR 1099 has higher surface fluxes than e Eri. Our model, 
which is displayed in Figure 13, is similar to Simon & Linsky’s 
(1980) model A, and should be accurate enough for our very 
schematic calculations. A full discussion of the validity of using 
a scaled € Eri model for HR 1099 is given by Simon & Linsky 
(1980). 

Initially, we assumed solar abundances (see, e.g., Morton 
1991) and a microturbulence model like that of Simon et al. 
(1980), in which the microturbulence rises from 1 km s _ 1 in the 
photosphere to 10 km s -1 for T > 8000 K. We used the radi- 
ative transfer code MULTI (Carlsson 1986) to solve for the 
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Fig. 13. — Model atmosphere for HR 1099’s K1 star based on the e Eri 
model computed by Simon et a!. (1980), and very similar to Simon & Linsky’s 
(1980) model A. The temperature (solid line, left axis) and electron density 
(dotted line , right axis ) are plotted as functions of mass column density. 


level populations of the H and Mg atomic models. For these 
schematic calculations we assumed complete redistribution 
(CRD) for hydrogen, although the method described by Carlss- 
on & Judge (1993) was used to mimic the effects of partial 
redistribution (FRD) on the hydrogen ionization. For the 
background opacities we included non-LTE bound-free hydro- 
gen, and all other sources were treated in LTE. For the radi- 
ative transfer calculations of the Mg n line profiles we adopted 
the PRD techniques of Uitenbroek (1989a, b). The Mg ii model 
used is that given by Harper (1992). The differences between 
the computed Mg ii resonance-line fluxes and those of Simon 
& Linsky (1980) are mostly due to the different treatment for 
hydrogen ionization. 

We computed line profiles for five different values of ft = 
cos 0, where 0 is the angle between the normal vector of the 
plane-parallel atmosphere and the observer’s line of sight. 
Three of the computed Mg ii k - line profiles are shown in 
Figure 14a. Based on the profiles computed for the five values 
of a disk-integrated profile is computed that allows for rota- 
tional broadening, assuming v sin i = 38 km s _1 (Vogt & 
Penrod 1983). After conversion from intensity to surface-flux 
units, the final disk-integrated Mg n k - line profile is compared 
with the observed profile in Figure 14 b. 

The integrated surface flux of the model k-line profile in 
Figure 14 b agrees well with the observations, and the model 
profile does have broad wings somewhat similar to those of the 
observed profile. However, the wings of the model profile are 
not broad enough. Also, the top of the model profile is too 
broad because the peaks of the self-reversed profiles in Figure 
14a are too far apart. Further computations, which retain the 
original electron density stratification, have revealed that 
lowering the microturbulence decreases the separation 
between the peaks of the disk-resolved self-reversed Mg ii pro- 
files. In Figure 15 we show the disk-integrated Mg ii h and k 
profiles that result when the microturbulence in the chromo- 
sphere is lowered from 10 to 3 km s _1 . The model fluxes in 
Figure 14 are a little too low, and this discrepancy becomes 
worse when the microturbulence is lowered, so we increased 
the magnesium abundance to 50% above solar in order to 



Fig. 14. — (a) Mg ii fc-line profiles of HR 1099’s K1 star computed from the 
model in Fig. 13 by use of the radiative transfer code MULTI (Carlsson 1986), 
assuming a solar abundance for magnesium and a microturbulence of 10 km 
s" 1 above T = 8000 K. Profiles are shown for three values of n = cos 9 , where 
9 is the angle between the normal vector of the plane-parallel model atmo- 
sphere and the line of sight. A disk-integrated profile is also shown (thick line), 
with corrections for rotational broadening. { b ) Disk-integrated profile from (a) 
compared with the observed profile. The interstellar absorption feature and the 
flux contribution of the G5 star (as determined by the multi-Gaussian fit in Fig. 
12) have been removed from the observed profile in ( b ). The velocity scale is in 
the rest frame of the K 1 star. 


produce model profiles in Figure 15 with fluxes that better 
match the observations. As in Figure 14 b, the bases of the 
theoretical profiles in Figure 15 are too narrow, so we have 
broadened the lines further by adding a macroturbulence of 40 
km s" 1 . The resulting profiles (dotted lines ) match the obser- 
vations better. 

Besides the h and k lines, we have also modeled the much 
weaker Mg ii subordinate line at 2798 A (centered at 270 km 
s _I in the top panel of Fig. 15). Our models do not show the 
subordinate line in emission when the opacities described 
above are used. Since this line is located in the wing of the 
Mg ii k line and we have not modeled the line overlap, we 
arbitrarily increased the contribution from the LTE opacities 
by a factor of 2 to mimic the increased opacity due to the k line, 
and the subordinate line then appeared in emission as shown 
in the top panel of Figure 15. Previous studies of this feature 
(e.g., Lemaire & Gouttebroze 1983) have shown that the for- 
mation of this line is expected to be sensitive to blending with 
the Mg ii resonance lines. This line has been observed in emis- 
sion above the solar limb (Feldman & Doschek 1977) and also 
in the spectrum of y Cru (M3.4 III), where fluorescence by Ly/? 
may be important (Carpenter, Robinson, & Judge 1995). 
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Fig. 15— Disk-integrated Mg n h ( bottom ) and k {top) line profiles com- 

puted as in Fig. 14, but with a microturbulence of 3 km s ' and with a 

magnesium abundance 50% higher than solar. The initial model profiles are 
shown as thick solid lines and are compared with the observed profiles that are 
plotted in histogram form. Interstellar absorption and the G5 star’s emission 
(as determined by the multi-Gaussian fits in Fig. 12) have been removed from 
the observed profiles. A better match between the computed and observed 
profiles can be obtained by adding a macroturbulence of 40 km s _l to the 
model profiles, which leads to the profiles shown as dotted lines. We have also 
modeled the Mg n subordinate line at 2798 A (centered at 270 km s 1 in the 
top panel), and the computed profile for this line is shown as a thin solid line. 
The fluxes between 200 and 300 km s' 1 have been increased by a factor of 10 
to better display the weak subordinate line. The velocity scale is in the rest 
frame of the K1 star. 


Since our model Mg n h and k line profiles have broad wings 
and fit the observed profiles fairly well when some macro- 
turbulence is added, it seems natural to propose that opacity 
and radiative transfer effects are the most important factors 
regulating the profiles of these lines. However, despite the simi- 
larity between the Mg n and C iv line profiles, it is very difficult 
to believe that similar processes could be controlling the C iv 
profiles. Solar C iv lines are limb-brightened (Kjeldseth-Moe & 
Nicolas 1977), implying optical depths much too low for 
opacity effects to induce broad emission wings. 

The optical depth of HR 1099’s C iv lines can be estimated 
by use of equation (5). In computing the observed surface flux 
F, we assume a stellar radius of 3.9 R 0 (Fekel 1983), and we 
assume that the C iv emission is uniform over the stellar 
surface. For Q = 6.4 (Cochrane & McWhirter 1983), g = 2, 
/= 0.1908 (Morton 1991), and log P e = 15.05 (see § 4.3), we 
find that x 0 — 0.72/AA for the stronger C iv line. 

It is important to be aware that AA in this equation is the line 
width associated only with the microturbulence. Only if 


opacity broadening and macroturbulence are unimportant can 
the observed width listed in Table 3 (FWHM = 194 km 
s“ 1 = 1.0 A) be used for AA. Otherwise, the observed width is 
only an upper limit for A A. A lower limit to the micro- 
turbulence present in the C iv line is that due to thermal 
motions, AA = 0.10 A. Using this lower limit, we calculate an 
upper limit for the optical depth of the C iv x 1 548 line of 
t 0 < 7.2. The only way that t g can be significantly larger than 
7.2 is if the C iv emission in reality originates from a very small 
area on the stellar surface, meaning we have grossly underesti- 
mated the true C iv surface flux. We believe that z 0 must be at 
least 1000 for opacity effects to produce broad wings, so it is 
highly improbable that the line wings of C iv are produced by 
radiative transfer effects. Considering also the very broad 
widths of the optically thin transition-region intersystem lines, 
the evidence for turbulence-dominated C iv profiles is over- 
whelming. For comparison, we note that equation (5) predicts 
r 0 ~ 10,000 for the Mg ii k line, which illustrates why our Mg n 
model profiles have opacity-broadened emission wings. 

Despite some success in reproducing the observed profiles of 
the Mg ii lines with radiative transfer models, we find it very 
unlikely that the C iv lines are also regulated by radiative 
transfer effects. Therefore, either the similarity between the 
Mg ii and C iv profiles is simply coincidental or the moderate 
success of the Mg n models is misleading and turbulence 
actually plays a much larger role in shaping the line profiles 
than the models suggest. To investigate this latter possibility 
further, we modeled the Si ii A1533 line, assuming a solar 
silicon abundance, PRD, the Si ii model used by Harper (1992), 
and a microturbulence of 3 km s” 1 (as in the Mg n model in 
Fig. 15). The result is shown in Figure 16. Our calculated 
surface flux for this line is too high by about a factor of 2, so we 
have divided the calculated fluxes by 2 in Figure 16 to allow for 



Fig. 16. — Observed Si u A 1533 line profile ( histogram ) and the single- 
Gaussian fit to the line and flat continuum assumed for the fit {thick solid lines). 
The velocity scale is centered on the centroid of the observed profile, which is 
redshifted relative to the K1 star by 13 km s 1 (see Table 3). The thin solid line 
is a Si n profile computed based on the atmospheric model in Fig. 13, assuming 
a solar silicon abundance and a microturbulence of 3 km s" '. The computed 
fluxes of this model were actually twice as high as those shown in the figure. 
We have divided the fluxes by 2 to allow for easier comparison of the com- 
puted and observed profiles. Adding a macroturbulence of 100 km s 1 to the 
model profile leads to the profile shown as a dotted line, which fits the data 
better. 
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Fig. 17.— Alternative line profile model to the two-component model in 
Fig. 9b. If the line profiles being produced by different regions of HR 1099’s 
atmosphere are Gaussian and have a normal distribution of widths (with a 
mean of 205 km s -1 and a = 100 km s' 1 ), the dotted line is the resulting 
disk-averaged line profile, which matches the observed profile {histogram) quite 
well. 


easier comparison between the modeled and observed Si n 
profiles. The model Si n line is found to have broad emission 
wings, but they are not as prominent or as broad as those of 
Mg ii. However, the observed Si n line profile is actually some- 
what broader than the observed Mg n lines, and the width of 
the model Si ii profile is therefore much too narrow. The 
observed line width can be reproduced only if a macro- 
turbulence of 100 km s ” 1 is added to the model profile, as 
shown in Figure 16. In such a model, macroturbulence actually 
plays a larger role in shaping the line profile than do the radi- 
ative transfer effects, suggesting that the K1 star’s chromo- 
sphere may in fact contain significant turbulence. Therefore, it 
at least seems possible that the broad wings of the Mg ii lines 
are in fact produced in part by turbulence. If the microflaring 
interpretation of the broad transition-region line wings is 
correct, then microflaring may also be occurring in the K1 
star’s chromosphere. We note that broad wings were not found 
in the chromospheric lines of Capella (Linsky et al. 1995b), so 
observations of other extremely active stars should be per- 
formed to determine how common broad chromospheric line 
wings are for active stars. 

4.6. Discussion 

Even if one accepts our assumption that the profiles of HR 
1099’s emission lines are regulated by turbulence, it remains 
open to debate whether microflaring is responsible for the par- 
ticularly broad line wings and whether the broad components 
of the two-Gaussian fits to these lines successfully quantify the 
microflaring contribution. It should come as no surprise that 
two-Gaussian fits, which imply two distinct emission com- 
ponents with different degrees of turbulence, are not unique 
models of these profiles. The flaring solar transition-region 
line profiles discussed by Cheng (1990) are very similar to the 
C iv profiles of HR 1099, which suggests that all of HR 1099’s 
transition-region line emission may be generated by micro- 
flaring. 

Figure 17 demonstrates that one can interpret HR 1099’s 
C iv profiles in terms of a single, normal distribution of line 


widths being emitted from the stellar surface. Assuming such a 
distribution may be reasonable because solar observations 
show that the widths of the C iv lines are not the same every- 
where on the solar disk (Dere & Mason 1993). Neither do all 
solar flares and transition-region explosive events produce 
identical line profiles (e.g., Cheng 1990; Dere et al. 1989). 

However, the solar example also suggests that a two- 
component model may be meaningful because the line profiles 
of transition-region explosive events are produced by a pheno- 
menon separate and distinct from the less dynamic processes 
that generate the narrower emission lines that are much more 
prevalent on the Sun. Therefore, it is not correct to claim that 
the widths of the very broad lines produced by the explosive 
events simply represent the high-velocity tail of a single solar 
transition-region line width distribution. It is more accurate to 
say that there are two separate distributions — one with rela- 
tively narrow line widths and one with broad widths. On the 
Sun, the narrow distribution completely dominates the disk- 
averaged line profile, but since the explosive events certainly 
appear to be a manifestation of stellar activity, this does not 
have to be the case on stars as active as HR 1099. It is natural 
to suppose that the explosive-event phenomenon may be a lot 
more prevalent on active stars, which nicely explains the broad 
line wings that have been observed in some of these stars, 
including HR 1099, and it explains why the lines of Procyon, a 
less active star, do not have the broad wings. However, we 
cannot rule out the possibility that differences in spectral type 
may produce effects that could also be at least partially 
responsible for the differences in the line profiles of HR 1099 
and Procyon. 

Besides the broad transition-region line wings, other 
observed properties of active stars are also suggestive of con- 
tinuous microflaring. On the Sun, coronal temperatures above 
10 7 K are observed only during solar flares, but very active 
stars such as HR 1099 have coronal temperatures this high 
even when they are not obviously flaring (see, e.g., Dempsey et 
al. 1993). The implication is that continuous small-scale flaring 
may be maintaining the very high coronal temperatures 
(Haisch, Strong, & Rodono 1991). Observations of active stars 
and solar flares have demonstrated that their X-ray lumi- 
nosities are strongly correlated with their microwave lumi- 
nosities. Because the same linear relation appears to apply to 
both the solar flares and the “quiescent” emission from the 
active stars, these studies suggest that the microwave and 
X-ray emission observed from active stars is produced by flare- 
like processes (Giide! 1994; Benz & Giidel 1994). 

All of this evidence suggests that microflaring is very preva- 
lent in the atmospheres of active stars. If the broad UV 
emission-line wings observed for HR 1099 and other active 
stars are in fact spectral signatures of this microflaring, then 
these UV line profiles give us a unique opportunity to 
study the properties of the microflares. Understanding this 
microflaring is very important, as many authors have proposed 
that microflaring may be the dominant source of coronal 
heating for the Sun and other stars (see, e.g., Porter et al. 1987; 
Parker 1988; references therein). 

5. SUMMARY 

In this paper we have analyzed high-quality GHRS obser- 
vations of both Procyon and HR 1099 in order to compare the 
atmospheric properties of a relatively inactive star with those 
of an extremely active star. The results of our analysis of the 
Procyon data are summarized as follows: 
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1. With the exception of the N v A1239 line, Procyon’s 
transition-region lines do not have the broad wings that have 
been observed on active stars, supporting our contention that 
the broad wings are probably a manifestation of stellar activity 
and may be caused by microflaring. Observations of other 
inactive stars are needed to determine conclusively whether 
inactive stars as a group do not have broad components in 
their transition-region lines. 

2. Procyon’s transition-region lines have redshifts very 
similar to the redshifts observed on the Sun. We find the red- 
shifts increase gradually with temperature, rising from +4 km 
s“ 1 2 at log T = 4.7 to -I- 1 1 km s' 1 at log T = 5.4. Line red- 
shifts such as these may indicate the presence of a circulation 
system in which most of the transition-region line flux is pro- 
duced by the downward leg of the system. 

3. Procyon’s strong transition-region lines all have blue- 
wing excesses. The cause of these excesses is uncertain, but 
perhaps they are spectral signatures of the upward leg of the 
aforementioned circulation system. 

4. The He n A 1640 line is probably formed in Procyon’s 
chromosphere since the line’s velocity and width are similar to 
those of other chromospheric lines but are not consistent with 
the transition-region lines. This suggests that the radiative 
recombination line formation model, rather than the col- 
lisional excitation model, is appropriate for Procyon. 

5. The O iv] density diagnostics indicate transition-region 
pressures of log P e < 15.5. This upper limit is inconsistent with 
the coronal pressures measured from EUVE data (log P e = 
16.0). The simplest explanation for this discrepancy is that 
Procyon’s atmosphere must be inhomogeneous. 

6. The emission lines of Procyon are significantly broader 
than those of the Sun. Except for O iv], Procyon’s transition- 
region lines all have essentially the same width, FWHM = 
78.7 ± 0.9 km s"\ which implies nonthermal velocities of 
£ = 46.1 ±0.5 km s 1 in the transition region. Procyon’s 
narrower chromospheric lines also have essentially the 
same width, FWHM = 39.6 ± 2.3 km s' 1 , which implies 
l = 22.7 ± 1 .6 km s " 1 in the chromosphere. 

7. The Si hi A 1207 line is the only transition-region line 
profile we have observed that is obviously altered by opacity 
effects. Our crude analysis of the Si iii profile suggests that 
t„ = 0.7 for the line, but we consider this to be only a lower 
limit to the actual optical depth. Based on this lower limit, we 
derive an upper limit for the electron pressure of log P e < 14.8. 
We note that this calculation relies on a stellar line surface flux 
computed under the assumption that the Si m emission is uni- 
formly distributed over the stellar surface. 

8. The Mg ii h and k line profiles of Procyon are very similar 
to those of the Sun. The line profiles are self-reversed and have 
stronger blue peaks, which are probably signatures of hydro- 
dynamical processes in the chromosphere. 

The results of our analysis of the GHRS observations of HR 
1099 are summarized as follows: 

1. The K.1 star accounts for nearly all of HR 1099’s UV line 
emission. Emission from the G5 star was detected only in the 
Mg ii h and k lines, where it is responsible for ~7.5% of HR 
1099’s Mg ii flux. It must be responsible for only ~5% or less 
of HR 1099’s C iv. Si iv, and He ii flux to explain the nonde- 
tection of the G5 star’s emission in those lines. 

2. Broad non-Gaussian wings were detected in four of HR 

1099’s strongest UV emission lines: Si iv A1394, C iv AA1548, 


1551, and He ii A1640. The broad components of the two- 
Gaussian fits to these lines account for 65.6% ± 1.0% of the 
total line flux. For comparison, the broad components of AU 
Mic and Capella account for ~40% of the line flux. Since HR 
1099 is more active than both AU Mic and Capella, the larger 
broad-component contribution for HR 1099 is consistent with 
a correlation between the broad components and stellar activ- 
ity and is consistent with the hypothesis that microflaring is 
responsible for the broad-component emission. After smooth- 
ing the data, we find that the N v A1239, A1 ii A 1 67 f , and Si in] 
A 1892 lines may also have profiles similar to those of Si iv, C iv, 
and He n. 

3. The broad components of the C iv and He ii lines are 
blueshifted relative to the narrow components. Linsky et al. 
(1995b) discovered that many of the transition-region lines of 
Capella have the same property. 

4. The transition-region lines of the K.1 star of HR 1099 
probably have intrinsic redshifts no larger than + 5 km s *, in 
contrast to the large 10-25 km s' 1 redshifts observed on the 
active G1 III star of Capella. 

5. The profiles of the K1 star’s C iv AA1548, 1551, Mg ii h 
and k , and He ii A1640 lines are very similar. Since the width of 
the He ii line is closer to that of the Mg ii lines than to the 
broader C iv lines, we have proposed that the He n line is 
formed in the chromosphere by radiative recombination. The 
extremely large X-ray surface flux of the K 1 star supports this 
conclusion. 

6. Using the density-sensitive O iv] lines, we derived a 
transition-region pressure (in units of cm 3 K) of log P e = 
15.05 ± 0.15. This pressure is similar to transition-region pres- 
sures measured for Capella and the quiet Sun. This is sur- 
prising, considering that the transition-region line surface 
fluxes of HR 1099 are about 4 times larger than those oi 
Capella and more than 100 times larger than those of the Sun. 

7. The S i and C i lines are much narrower (FWHM « 5C 
km s' l ) than the widths of lines formed in the upper chromo- 
sphere and transition region (FWHM = 173 ± 19 km s *). 
which suggests a large increase in turbulence from the lower to 
the upper chromosphere. 

8. We see no obvious evidence of opacity in any of the line 
profiles of HR 1099. The extremely broad widths of the opti- 
cally thin intersystem lines suggest that nonthermal velocities 
control the transition-region line profiles. 

9. The profiles of the C iv and Mg n lines are remarkabh 
similar, which suggests that similar processes regulate the pro- 
files of these lines. It is possible that microflaring may be 
responsible for the broad wings observed in both the C iv and 
Mg ii lines. This would imply that, on HR 1099’s K1 star 
microflaring may be taking place in the chromosphere as wel 
as in the transition region. However, radiative transfer modeF 
suggest that, for Mg ii, opacity effects can in principle also 
account for the broad wings. 
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